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Foreword 


This  is  part  I  of  the  final  technical  report  for  Contract  No.  N62269- 
76-C-0378,  which  is  sponsored  by  the  Naval  Air  Development  Center,  Warminster, 
Pa.  The  work  was  performed  during  the  period  of  July  1,  1976  through 
December  30,  1977.  Mr.  Lee  W.  Cause  was  the  contract  monitor. 

The  contracted  study  is  iinder  the  title  "Certification  of  Composite 
Aircraft  Structures  under  Impact,  Fatigue  and  Environmental  Conditions"; 
parts  I  and  II  of  the  study  are  under  the  supervision  of  Dr.  P.C.  Chou,  while 

III  is  xmder  Dr.  A.S.D.  Wang,  both  of  Drexel  University. 

This  report  concerns  part  I  of  the  contract,  low  speed  impact  of  plates 
of  composite  materials.  It  is  a  self-contained  report,  including  definitions 
of  all  nomenclature  used,  and  its  own  introduction  and  conclusions. 

The  authors  would  like  to  thank  Dr.  Edward  J.  McQuillen,  Dr.  James  L. 
Huang  and  Mr.  Lee  W.  Cause  for  the  frequent  technical  discussions.  They 
would  also  like  to  thank  Mr.  Frank  Patota  and  Mr.  Ceorge  Chou  who  helped 
conducting  part  of  the  experiments. 
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I.  INTRODUCTION 

Aerospace  structures  frequently  undergo  impacts  by  blunt  objects,  in-- 
eluding  dropped  tools,  hail,  and  runway  stones  impinging  on  exposed  aircraft 
components,  and  foreign  objects  entering  jet  engines.  The  failure  of  a 
component  subject  to  such  an  impact  is  often  of  the  structural  type,  rather 
than  due  to  local  penetration  or  indentation.  A  number  of  numerical  techniques 
(e.g.,  f inite-^element  method,  lumped-parameter  models)  are  available  for  calcu¬ 
lating  structural  response  to  impact,  but  these  are  typically  complicated, 
time-consuming,  and  usable  only  on  a  problem-by-problem  basis.  Designers  need 
a  quick,  convenient,  widely  applicable  method  for  this  purpose. 

In  a  previous  report  [1],  a  method  was  developed  for  constructing  a  de¬ 
sign  curve  which  predicts  the  response  of  a  given  type  of  structure  to : impact 
loading.  This  curve  gives  the  maximum  strain  in  the  structure,  which  may  have 
various  dimensions  and  material  properties,  due  to  impacts  involving  a  certain 
range  of  impact  masses  and  velocities.  An  example  of  the  bending  response  of 
a  simply  supported  beam  under  central  impact  was  presented  in  detail.  Both 
experimental  results  and  numerical  calculations  involving  several  solutions 
of  beam  impact  were  used  in  establishing  the  design  curve.  The  impact  cases 
studied  in  [1]  were  limited  to  large  impactor  mass,  where  the  impactor  has 
mass  roughly  equal  to  or  greater  than  that  of  the  beam. 

In  the  present  report,  the  design-curve  approach  is  extended  to  anisotropic 
plates,  and  impacts  by  small  impactors.  Shear  failures  are  also  studied.  The 
elementary  model  of  impact  on  beams  presented  in  [1]  is  generalized  to  embrace 
all  structures.  ^This  model  is  then  applied  to  the  cases  of  clamped  and 
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simply  supported  anisotropic  rectangular  plates  and  to  predicting  impact 
failure  due  to  shear  effects  in  both  beams  and  plates.  Design  curves  are 
also  developed  for  predicting  the  response  of  beams  and  plates  to  impact 
by  small  impactors.  In  each  Instance,  both  analytical  methods  and  Impact 
experiments  are  employed  in  generating  the  design  curves. 

In  dealing  with  low-velocity  impact  problems,  it  is  convenient  to  divide 
them  into  two  domains,  based  on  the  ratio  M  of  the  mass  of  the  structure  to 
the  mass  of  the  impactor  (see  Fig.  1).  The  motion  of  the  impactor  and  the 
response  of  the  structure  for  a  large  impactor  differ  greatly  from  those  for  a 
small  impactor.  The  two  domains  will  be  treated  separately;  the  assumptions 
and  final  design  curves  are  also  different  for  the  two  domains. 

In  Impacts  where  the  mass  ratio  M  is  small  (large  impactor) ,  it  has 
been  observed  that  multiple  collisions  occur  and  that  the  general  motion  of 
the  Impactor  follows  the  path  of  the  structure  at  the  impact  point  (Fig.  la). 
The  final  rebound  of  the  impactor  takes  place  after  the  structure  has  reached 
its  maximum  deflection.  Further,  the  entire  event  is  generally  more  pro¬ 
longed  than  the  fundamental  period  of  structural  vibration. 

On  the  other  hand,  when  M  is  large  (small  impactor),  only  a  single, 
sudden  collision  occurs,  after  which  the  impactor  rebounds  and  the  structure 
continues  to  vibrate  freely,  reaching  its  maximum  deflection  at  a  later  time 
(Fig.  lb).  The  duration  of  actual  contact  is  characteristically  much  shorter 
than  the  fundamental  period  of  vibration.  The  sudden  rebound  of  the  Impactor 
is  due  chiefly  to  the  elastic  resistance  of  the  structure  and  the  Impactor 
to  local  indentation,  i.e.,  contact  effects  between  the  two  bodies. 
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In  [1],  six  analytical  models  of  impact  were  examined,  including  two 
classical  one-degree-of-freedom  models;  a  two-degrees-of-freedom  model  which 
accounts  for  contact  effects;  a  solution  by  Clebsch  [2],  which  assumes  plastic 
impact  and  treats  the  beam  with  attached  impactor  as  a  free-vibration  problem; 
a  modified  Clebsch  solution  due  to  McQuillen  et  al.  [3];  and  Timoshenko's 
solution  [4],  which  couples  the  impactor  displacement  with  the  beam  deflection 
using  Hertz's  law  for  contact  deformations.  The  dimensionless  parameters 
which  determine  the  peak  impact  response  according  to  each  solution  method 
were  derived. 

It  was  shown  in  [1]  that,  for  those  models  which  neglect  contact  effects, 
the  maximum  strain  due  to  impact,  when  properly  normalized  to  £,  is  a  function 
of  the  mass  ratio  M  only  and  is  independent  of  any  other  parameters.  These 
solution  methods,  however,  are  generally  not  suited  for  treating  problems 
involving  small  impactors.  Conversely,  the  two  models  that  include  contact 
effects  indicate  a  dependence  of  e  on  M  and  another  parameter,  and  are  useful 
over  the  entire  range  of  M.  For  large  impactors,  it  was  demonstrated  that  the 
dependence  on  the  other  parameter  is  weak,  so  that  a  single  e  vs.  M  curve 
approximately  represents  all  impact  situations.  In  addition,  when  data  from 
assorted  impact  experiments  were  plotted  in  coordinates  of  e  and  M,  a  single 
curve  could  be  drawn  through  all  data  points,  with  about  +  30%  variation,  for 
large  impactors. 
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In  this  report,  several  approaches  are  taken  in  extending  the  design- 
curve  concept  to  impacts  of  other  structures  by  large  impactors  (Section  II). 
First,  the  one-degree-of-freedom  model  described  in  [1]  is  generalized  to 
encompass  all  structures  and  is  then  applied  to  anisotropic  laminated  rectang¬ 
ular  plates  with  clamped  and  simply  supported  boundary  conditions.  Then,  in 
order  to  identify  additional  significant  parameters,  the  structural  impact 
equations  of  Timoshenko  are  cast  in  a  dimensionless  forny  the  influence  of 
each  of  these  parameters  is  appraised  by  a  series  of  calculations  using 
Timoshenko's  solution  in  which  the  values  of  the  parameters  are  systematically 
varied.  Next,  we  discuss  two  series  of  impact  experiments  performed  on  lam¬ 
inated  graphite /epoxy  plates;  in  one  series,  the  bending  strains  in  the  plates 
were  measured,  and  in  the  other,  the  plates  were  repeatedly  impacted  at  grad¬ 
ually  increasing  velocities  until  failure  occurred.  The  purpose  of  these  ex¬ 
periments  is  twofold:  to  verify  the  analytically  constructed  design  curves,  and 
to  demonstrate  how  design  curves  may  also  be  constructed  experimentally.  Finally, 
design  curves  for  predicting  shear  failure  due  to  impact  are  developed,  and 
the  relative  importance  of  shear  and  bending  effects  in  impact  failure  is  dis¬ 
cussed. 

For  cases  in  which  the  mass  of  the  Impactor  is  small  with  respect  to  the 
mass  of  the  structure  (large  M) ,  the  dependence  of  the  generalized  strain  e 
on  the  mass  ratio  M  is  not  exclusive;  another  dimensionless  parameter  which 
involves  the  geometric  and  material  properties  governing  the  contact  effects 
between  the  two  bodies  can  also  be  shown  to  be  significant.  This  contact 
parameter,  denoted  X,  and  the  mass  ratio  M  together  determine  the  duration  of 
contact  between  the  structure  and  impactor  (relative  to  the  fundamental  period 
of  vibration),  which  is  typically  quite  short.  Clearly,  the  duration  of  contact 
is  important  in  determining  the  response  of  the  structure. 
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The  most  sophisticated  solution  method  applicable  to  problems  involving 
small  impactors,  the  Timoshenko  solution,  involves  a  nonlinear  integral 
equation  which  can  only  be  solved  numerically  on  a  problem-by-problem  basis, 
and  is  thus  not  convenient  for  developing  a  design  curve.  However,  by  making 
the  simplifying  assumption  that  the  structure  does  not  appreciably  deflect 
the  short  period  of  actual  contact,  the  problem  uncouples  into  two 
parts  —  the  elastic  Hertzian  impact  of  a  sphere  on  a  semi— infinite  body,  and 
the  vibration  response  of  the  structure  to  a  dynamic  load  (i.e.,  the  contact 
force) .  This  leads  to  an  approximate  relationship  (for  each  type  of  structure) 
between  the  parameters  e,  M,  and  X,  which  may  be  further  simplified  to  a 
simple  equation  in  the  mass  ratio  M  and  a  new  generalized  strain  for  small 

/ 

impactors  e*  (Section  III) .  By  comparison  with  experimental  data  and  with 
calculations  involving  the  complete  Timoshenko  solution,  the  form  of  this  re¬ 
lationship  is  shovm  to  be  useful  as  a  design  guide. 

In  Section  IV,  procedures  for  using  the  design  curves  presented  in  this 
report  are  described,  and  a  few  examples  are  given.  Also  included  are  methods 
for  constructing  new  design  curves  for  structures  not  treated  in  this  report, 
based  on  either  analytical  tools  or  impact  experiments. 
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II.  DESIGN  CURVE  FOR  IMPACTS  BY  LARGE  IMPACTORS 

For  impact  cases  where  the  mass  of  the  structure  is  small  compared  to 
the  mass  of  the  impactor,  the  structure  may  be  assumed  to  behave  as  a  simple 
spring,  or  as  a  one-degree-of-freedom  mass-spring  system.  Details  of  this 
approach  as  applied  to  simply  supported  beams  were  presented  in  Ref. [1];  in 
this  section,  this  impact  model  is  briefly  reviewed,  becatise  it  will  be  useful 
in  the  development  of  parameters  in  plate  impact  problems.  These  parameters 
are  employed  in  constructing  a  design  curve  for  predicting  the  response  of 
plates  to  impacts  by  large  masses.  Finally,  by  examining  data  from  a  large 
number  of  experiments  and  by  studying  impact  solutions  using  more  sophisticated 
analytical  techniques,  the  validity  of  the  design  curve  is  established  and  the 
range  of  its  usefulness  is  demonstrated. 

A.  Generalized  One-Degree-of-Freedom  Model 

In  this  model,  the  impactor  and  the  structure  are  considered  attached  to¬ 
gether  after  contact  and  move  together  as  a  single  mass  m,  and  the  motion  of 
the  structure  is  governed  by  an  equivalent  spring.  The  initial  velocity  of 
the  combined  mass,  Vq,  may  be  determined  by  two  methods:  One  is  based  on 
the  conservation  of  momenttna  (Fig.  2),  the  other  on  the  conservation  of  energy 
(Fig.  3).  In  the  first. case. 

m  =  m2  +  em^^  (1) 

where  m2  is  the  impactor  mass,  and  em^^  is  an  "equivalent"  mass  of  the  structure. 
The  equivalent  mass  may  be  obtained  by  matching  the  total  kinetic  energy  of  the 
structure  with  the  kinetic  energy  of  the  unknown  equivalent  mass  traveling  at  the 
velocity  of  the  impact  point  in  the  structure,  assuming  that  the  deflection  mode 
shape  is  the  static  deflection  cvO've.  For  simply  supported  beams, 
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(a)  beiora  impact  (b)  Iramcdiatcly  after  impact 

Fir’*T^c  2m  _  One-degree<-of-‘freedom  model  of 
impact,  momentum  conserved# 


(a)  before  impact  (b)  at  maximum  deflection 

Figure  :3«  Onc-degrcc-of-freetiom  model  of 

impact,  energy  conserved* 
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e  =  17/35.  By  equating  the  momentum  of  the  impactor  before  impact,  m2  v, 
with  the  momentum  of  the  combined  mass  m,  we  have 


where 


Vq  =  v/ (1  +  eM) . 
M  = 


Note  that  according  to  Eq.  (2),  the  energy  before  and  after  impact  is  not 

2  2 

conserved,  i.e.,  mv^  • 

In  the  energy  conserved  case,  we  assume  e  =  0,  so  that 


Vq  =  V  (4) 

As  a  result,  the  kinetic  energy  in  the  impactor  is  conserved,  and  will  be 
entirely  converted  to  strain  energy  in  the  structure.  For  small  values  of 
M,  the  difference  between  these  two  approaches  is  small. 

The  spring  constant  is  the  force  per  unit  deflection,  with  respect 
to  a  force  acting  at  the  impact  point  in  the  direction  of  the  impactor 


velocity,  or 


P  =  Ki  Wi 


After  acquiring  the  initial  velocity  Vq,  the  mass-spring  system  is  assumed 
to  perform  a  free  vibration.  The  maximum  deflection  is  then 

In  modeling  the  structure  by  a  mass-spring  system,  it  is  implied  that  only 
the  first  mode  of  vibration  is  retained  and  that  its  mode  shape  is  the  same 
as  the  static  deflection  distribution  tinder  a  concentrated  force. 


Next,  we  shall  consider  the  most  critical  strain  the  structure  which 

occurs  at  a  known  point  2.  Assume  that  a  proportional  relation  between  £2  and 
Wj^  can  be  found. 
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where  dj^2  determined  from  the  static  deflection  distribution,  or  from 
experimental  static  measurements.  Combining  Eqs.  (5)  and  (6),  we  obtain  the 
maximum  strain  as 

^0  1  /  m 


=  0  1  m 

■  112  K 


or,  combining  Eqs.  (1),  (2),  and  (7), 


'2  “  d,„ 

max  12 


(l+eM)K, 


1.  Simply  Supported  Beam 


The  next  step  of  defining  a  generalized  strain  e  is  best  illustrated  by 
considering  a  specific  type  of  structure.  For  a  simply  supported  beam  of 
length  L  and  depth  h.  Impacted  at  middle  span,  we  have 

d^2  =  L^/6h  (9) 


=  48EI/L^  . 


Therefore, 


e  =  JLlL\  /IZ^l 

max  k  y4  M(l+eM) 

where  c^  =  /e/p,  the  velocity  of  longitudinal  waves  in  a  bar. 
Defining  the  generalized  strain  as 


e  =  E 


we  obtain  finally. 


max  V  h 


^  1/4  M(l+eM) 


which  is  the  equation  of  the  design  curve  for  simply  supported  beams  impacted 
by  large  masses,  as  given  in  [1]. 
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I 


I. 

1 


The  corresponding  expressions  for  ^  simply  supported 

and  clamped  orthotropic  plates  will  be  given  in  the  following  sections. 


2.  Simply  Supported  Plate 

The  impact  design  curve  for  simply  supported,  rectangular, 
orthotropic  plates  may  be  generated  by  applying  the  generalized  one-degree- 
of-freedom  model.  For  transverse  impacts  at  the  center  of  the  plate,  we  first 
consider  the  corresponding  static  problem  Involving  a  centrally  applied  con¬ 
centrated  load.  In  this  case,  the  maximum  deflection  occurs  at  the  center  of 
the  plate,  and  the  maximum  strain  on  the  surface  opposite  the  applied  load. 

It  is  shown  in  Appendix  A  that  subject  to  a  mild  approximation,  the  spring 
constant  of  a  plate  so  loaded  may  be  expressed  as 


4  _ 
IT  D 


11 


^  '^max  Aab  f^(n) 


where 


I 


f-.(n)  ®  i,  ^  c 

^  m=l,3,5  n=l,3,5  W 


00 

I  7^ 


(14) 


(15) 


C  =  m^  +  2m^n^n  +  n^ri^ 
mn 


(16) 

(17) 


The  strain-deflection  constant  relating  the  strain  due  to  bending  in 
the  x-direction  to  the  maximum  deflection  Wj^  is 


e  T  TT^h  f„(Ti) 

a  -L.  a  2 

''l  ^12  2ab  (f)fi(n) 


(18) 


11 
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fo(n)  =  I 


m“l,3,5  n=l,3,5  mn 


By  combining  Eqs.  (14)  and  (18)  with  Eq  *  (8),  and  assuming  that  e  =  0 
(conservation  of  energy),  we  obtain 


faCn)  hv 


hv__v:iL_ 

[M 


If  we  define  the  generalized  strain  e  as 

X 


e  =  e 

X  X 


where  Cp  -  the  speed  of  flexural  waves  in  the  plate  in  the 

x-direction,  and  k  is  the  radius  of  gyration  of  the  plate,  then, 

*  iTm" 

where 

gj^(n)  =  f2(n)/7f^Cn) 

This  function  is  plotted  in  Fig.,  4.  Similarly,  we  have 

g,  (h’b 


Equations  (22)  and  (24)  may  be  plotted  as  Impact  design -cm: ves ,  as  shown 
in  Fig.  6.  Note  that  in  this  case,  the  generalized  strains  e  and  e 

X  y 

are  functions  not  only  of  the  structure-to-impactor  mass  ratio  M  but  also 
of  the  aspect-orthotropy  ratio  q . 


3.  Clamped  Plate 

We  may  develop  the  design  curve  for  clamped-edge ,  rectangular, 
orthotropic  plates  in  a  similar  manner.  The  spring  constant  of  such  a  plate 

i 

with  respect  to  a  central  point  load  may  be  expressed  approximately  (see 


12 


NADC-7fc259  -60 


Appendix  B)  as 

bD 

K  =  •  ■  ■  (25) 

^  a-^  £3(0) 

The  strain-deflection  constant  d^2  relating  the  strain  in  the  x-direction  with 


the  maximum  deflection  is 


a2  £3(0) 

^12  "  h  f^(Ti) 


where  £3(0)  and  f^Cn)  are  given  in  Appendix  B. 

By  combining  Eqs.  (25)  and  (26)  with  Eq.  (8),  and  assuming  that  e  =  0 
(conservation  of  energy),  and  defining  the  generalized  strain  according  to  Eq.  (21) 


we  obtain 


82(11) 


where 


g,(n)  =  f^(n)//f '(n) 


This  function  is  also  plotted  in  Fig.  4.  Similarly,  we  have 


Equations  (27)  and  (29)  represent  the  impact  design  curves  for  a  clamped 
rectangular  plate.  Again,  the  generalized  strains  e  and  e  are  functions 

X  y 

of  both  the  mass  ratio  M  and  the  aspect-orthotropy  ratio  ri. 

4.  Comprehensive  Impact  Strain  Curve 

Examination  of  the  design-curve  equations  developed  above,  Eqs.  (22), 
(24),  (27),  and  (29),  shows  that  a  different  vs.  M  curve  must  be  prepared 
not  only  for  each  boundary  condition  but  also  for  each  value  of  the  aspect- 
orthotropy  parameter  p.  These  curves  can  be  combined  into  a  single,  com¬ 
prehensive  design  curve  if  we  introduce  a  new  dimensionless  generalized 
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Figure  4.  Functions  of  aspect-orthotropy  parametex-  rtj  for  constructing  inqpact  design 
curves  for  clai!5>ed  and  simply  supported  rectangular  orthotropic  plates. 
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Strain,  e,  defined  as 


e  = 


X 


g,.  (n) 


,  i  =  1,2 


(30) 


where  applies  to  simply  supported  plates  and  is  given  by  Eq.  (23)  and 


§2(0)  applies  to  clamped  plates  and  is  given  by  Eq.  (28).  Thus,  the 


equation  for  the  design  curve  in  e  vs.  M  coordinates  becomes 

7  =  (31) 

All  impact  data,  regardless  of  the  boundary  condition  (simply  supported  or  clamped) 
or  the  value  of  the  parameter  n  may  be  presented  on  a  single  curve.  This 
curve  is  shown  in  Fig.  5,  in  which  are  included  data  from  numerous  impact 
experiments  which  are  discussed  in  detail  in  a  later  section. 


5.  Critical  Impact  Velocity  Curve 

The  design  curves  in  terms  of  or  e  vs  M  give  the  maximum  strain  in  the 
plate  for  a  given  impact  situation.  This  strain  may  be  compared  to  the 
ultimate  failure  strain  of  the  plate  (assuming  a  maximum- strain  failure 
criterion)  to  determine  whether  the  plate  fails  as  a  result  of  a  given  impact. 
However,  by  introducing  a  new  parameter,  the  dimensionless  impact  velocity, 
we  can  construct  an  alternate  form  of  the  design  curve  which  can  be  directly 
used  to  determine  whether  a  given  impact  causes  failure  of  the  plate.  Again, 
the  curve  may  be  developed  either  by  analysis  or  by  experiment.  To  do  this. 


a  dimensionless  impact  velocity,  v,  is  defined  as 

vh 


V  = 


(32) 


where  v  is  the  impact  velocity,  h  is  the  plate  thickness,  Cp  is  the  speed  of 


flexural  waves  in  the  plate,  and  k  is  the  radius  of  gyration  of  the  plate. 
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Figure  5.  Comprehensive  impact  strain  curve  for  simply  supported  and 
clamped  rectangular  plates. 
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« 


1 


Then,  by  Eq's  (21),  (22)  and  (32)  the  critical  dimensionless  velocity  v^. 

(i.e.,  the  lowest  value  of  v  at  which  failure  occurs)  is  related  to  the  failure 
strain  by 


'^c  “  ^f  g.(n) 


(33) 


In  V  vs.  M  coordinates,  the  curve  v  =  v  (M)  divides  the  plane  into  a  safe 

c 

region  (below  the  curve)  in  which  failures  due  to  Impact  do  not  occur  (according 
to  the  assumed  maximum-strain  failure  criterion)  and  an  unsafe  region  (above 
the  curve)  in  which  failures  do  occur.  That  is,  if  the  point  (v,M)  corresponding 
to  a  given  impact  situation  falls  above  the  curve,  then  failure  occurs;  if 
the  point  falls  below,  no  failure  occurs.  Thus,  Eq.  (33)  is  useful  as  a  de¬ 
sign  curve  equation.  A  typical  v^  vs.  M  curve  is  shown  in  Fig.  12. 

Further,  if  we  define  a  new  dimensionless  impact  velocity  v  such  that 

V  =  g^(h)  V 

then  the  equation  of  the  alternate  design  curve,  in  terms  of  v  vs.  M,  is 


V  =  E-  (34) 

c  f 

This  single  curve,  v  =  v^(M),  may  be  used  to  predict  impact-induced  failure 
for  all  plate  structures,  regardless  of  the  boundary  conditions  or  the  value 
of  the  aspect-orthotropy  ratio  n  (provided  these  structures  have  the  same 
value  of  failure  strain  e^) .  Note  that  the  value  of  can  be  taken  as  the 
nominal  value  (say,  the  static  ultimate  strain)  of  the  given  material  or  can 
be  determined  by  performing  a  single  impaCt-to-f allure  test.  A  typical  v^ 
vs.  M  curve  is  shown  in  Fig.  13. 

B.  Timoshenko  Solution  of  Transverse  Plate  Impact 

Timoshenko's  approach  for  solving  transverse  impact  problems  on  simple 
beams  by  coupling  Hertz's  law  of  contact  with  the  Euler  beam  equation  has  been 
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extended  to  the  case  of  a  simply  supported  rectangular  plate  by  Eringen  [5] 
and  more  recently  to  the  case  of  a  simply  supported  anisotropic  laminated 
plate  by  Sun  and  Chattopadhyay  [6].  In  this  section  we  review  the  salient 
points  of  the  latter  solution  and  show  how  it  may  be  used  to  develop  an  impact 
design  curve  for .anisotropic  plates.  In  addition,  a  normalized  form  of  the 
equations  corresponding  to  a  special  case  of  this  solution  is  derived,  so 
that  the  dimensionless  parameters  governing  the  peak  strain  response  may  be 
identified.  Finally,  a  series  of  numerical  computations  is  performed,  in 
which  each  of  the  derived  parameters  is  systematically  varied  in  order  to 
demonstrate  the  significance  of  each  parameter  in  describing  the  response  to 
a  given  Impact. 

1.  Description  of  Solution  Method 

The  anisotropic  plate  equations  of  Whitney  and  Pagano  [7]  are  used 
to  predict  the  plate  motion.  Neglecting  effects  of  rotatory  inertia,  the 
deflection  of  a  symmetric  cross -ply  laminated  plate  due  to  a  centrally  applied 
force  F(t)  is 

nrhi-Z 

4  f(-l)  ^  f*" 

w(x,y,t)  =  —  EE  — - -  F(T)sina)  (t-T)dT 

“l  “  n  I  Jq  mn 


/ 
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Where  is  a  constant  and  a  Is  the  indentation  of  the  impactor  relative 
to  the  plate  surface,  or 

^  “  =  W2  -  w  (|,  I",  t) 

for  central  impact,  where  W2  is  the  Impactor  displacement. 

Newton's  law  applied  to  the  impactor  m2  is 


w-  =  vt  - 

I  m. 


‘2  JQ  Jo 


F  dt  dt 


These  four  equations  (35-38)  may  be  combined  into  a  single  nonlinear  integral 
equation  in  terms  of  the  contact  force  F  between  the  plate  and  the  Impactor, 


=  vt  - 


‘2  JQ  JQ 


F  dt  dt 


/ 

4  ^  “  ft 

~  I  I  ~~  FCt)  sin  to  (t-T)di 

“1  m=l,3,5  n=l,3,5  “mn  Jq  ^ 


which  is  solved  numerically  by  applying  the  small-increment  method  suggested 
by  Timoshenko,  in  which  the  contact  force  F  is  assumed  to  be  constant  or 
linearly  varying  during  any  time  increment  At.  Expanding  the  above  integrals 
to  calculate  the  force  F  during  the  i*^^  time  interval,  we  obtain 


=  v,At - SD..,,F. - F 

i  “2  j=l  J  “1  j=l  j 


where 


m  n  0) 


■j  {cos[a)^^(i-j)AT] 


If  the  contact  force  is  approximated  as  a  piecewise  linear  continuous 
function  of  time,  with  an  average  value  of  F^  during  the  1^^  time  step,  then 
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ji  “i-j+l  ^3  ■  (i-2>(F2-Fi) 

+  (i-3)  (F3-F2+F^)+“*+(i-j)(F^-Fj_^+Fj_2-***+  F^)+‘ •  •+ 


For  computing  the  solution  of  Eq.  (40),  Sun  and  Chattopadhyay  have 
suggested  a  recursion  method,  but  we  have  found  that  such  a  time-saving 
^PP^ox^^^tion  is  unnecessary .  A  listing  of  the  computer  program  for  solving 
this  equation  is  presented  in  Appendix  C. 


2 •  Dimensionless  Form  of  Timoshenko.  Solutions 
a,  Plate  Impact 

By  expressing  the  Timoshenko  solution  in  a  particular  normalized 
(dimensionless)  form,  the  parameters  governing  the  impact  response  of  plates 
may  be  identified.  Specifically,  it  is  shown  that  the  generalized  strain  e 

X 

depends  on  three  dimensionless  parameters:  the  mass  ratio  M,  the  aspect- 
orthotropy  ratio  n,  and  also  another  parameter  which  involves  the  geometric 
and  material  properties  governing  the  contact  effects  between  the  plate  and 
the  impact or. 

For  a  simply  supported  rectangular  plate  impacted  at  the  center, 
let  us  consider  the  governing  equation  (39)  and  the  strain  equation. 


; 
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_  2hir  ^ 
~  2  ^ 


2  ft 
Z  —  f 


F(t)  sin  0)  (t-T)dT 
nxi 


If  we  assume  that  the  plate  is  specially  orthotropic  in  bending  (i.e., 
=  0)  and  further  that  the  flexural  rigidities  are  related  by 


^12  + 


(see  Appendix  F)  then  the  natural  frequencies  are 


where 


w  “  “n 
mn  11  mn 


Uh  ; 


(1+n) 


A  (n)  = 
mn 


2  „  2 
m  +  nn 


Note  that  this  expression  for  the  natural  frequencies  is  a  special  case 
of  that  employed  in  the  Sun  and  Chattopadhyay  solution,  in  which  the 
approximation  of  Eq.  (42)  is  not  used. 

In  order  to  normalize  Eqs.  (39)  and  (41),  the  following  dimensionless 


variables  are  defined. 


t  . 

P  -  F 


Recalling  the  definition  of  generalized  strain, 

-  .  -iSiiax  ^  y. 

X  hv  P 

then  Eqs.  (39)  and  (41)  may  be  written  as 


X  1/3  p2/3  ^  -  _  .1 

P  /(1-fn)^ 


F  dt  dt 


— T  ^  ^  a"^'/  \  f  F(t)  sin  A  (t-T)dT 

/  m  n  A^(n)  L  ^  mn 
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2  ft 


“  ~Z - T  ^  ^  A — F(T)sin  A  (t-T)dT 

TT^i+n)^  Jo  ^ 


where 


M  = 


“ll'^V  “u 


and  T  is  a  dummy  variable  for  dimensionless  time  t.  It  can  be  seen  from 

these  normalized  equations  that  the  generalized  strain  depends  only  on 

the  mass  ratio  M,  aspect-orthotropy  ratio  n,  and  the  parameter  A  ;  that  is 

P 

The  parameter  defined  in  Eq.  (48)  may  also  be  expressed  as 


Thus,  the  parameter  A^  can  be  described  as  the  product  of  several  dimension¬ 
less  quantities. 

^  ^  /  elastic  force\^  wave  speed 

p  an  )  X  contact  force  /  ^  \  impact  velocity /’ 

X  t e ^leng th^ ^  [function  of  aspect-orthotropy  ratio] 


b.  Beam  Impact 

In  Ref. [l],  the  equations  for  the  Timoshenko  solution  of  beam  impact 
were  presented  and  discussed.  Tliese  equations  may  also  be  normalized  by  the 
above  procedure.  If  we  define 


t  =  co^t 
F  =  FL^oj^/EIv 
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then  the  equations  governing  the  response  of  a  beam,  corresponding  to 
Eqs.  (39)  and  (41)  for  the  plate,  may  be  written  as 


f2/3  .  - 


-JL  r 
/  io 


ft  _  _ 

F  dt  dt 

•'o 


and 


where 


A. 

.^1  .2 


F(t)  sin  i^(t-T)dT 


e  =  ^  E 
1  } 


F(t)  sin  i^(t-T)dT 


(51) 


(52) 


=  /Eiy  ^  ^  /  k\ 


which  may  be  expressed  as 


2  /  Kj  h 
2^ 


=  (constant  x 


^deflection  forc^^ 


contact  force 


X 


wave  speed 


radius  of  gyration 
depth  of  beam 


ion^ 


(53) 


^impact  velocity^ 

It  can  be  seen  from  Eqs.  (51)  and  (52)  that  in  this  case  the  generalized 

strain  e  depends  only  on  the  mass  ratio  M  and  the  contact  parameter  X.  ,  or 

b 


e  =  e(M,Xj^) 


In  Ref.  [1]  it  was  shown  that  for  large  impact  masses  (M  <  ~  2)  the 

dependence  of  the  generalized  strain  e  on  the  contact  parameter  X,  was  only 

b 

slight;  that  is 


e  =  £(M)  for  small  M. 

For  the  case  of  plate  Impact,  the  influence  of  the  contact  parameter  X^  is 
also  relatively  insignificant  in  determining  the  peak  strain  response  to 
impacts  by  large  masses.  This  will  be  demonstrated  in  the  following  section. 
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3,  Parametric  Study  for  Large  Impactors 

As  can  be  seen  from  Eq.  (49),  the  response  of  a  specially  orthotropic 
plate  satisfying  the  condition  of  Eq,  (42)  is  described  by  several  param¬ 
eters:  M,  n,  and  In  this  section,  we  examine  the  significance  of  each 

of  these  parameters  for  cases  where  the  mass  of  the  impactor  is  large  com¬ 
pared  to  the  mass  of  the  plate.  It  will  be  demonstrated  that,  for  plates 
of  given  aspect-orthotropy  ratio,  n,  the  dependence  of  the  impact  response 
on  parameters  other  than  M  is  weak  for  low  values  of  M,  so  that 

e  ^  e  (M,n)  for  small  M, 

X  X 

To  determine  the  significance  of  the  various  parameters  to  the  impact 
response,  several  plate  impact  problems  have  been  calculated  using  the 
Timoshenko  solution  method.  Each  problem  of  this  series  is  based  on  an 
impact  situation  with  the  impact  constants  listed  in  Appendix  C,  A  para¬ 
metric  study  was  carried  out  by  systematically  varying  the  values  of  the 
impact  velocity  v,  the  Hertzian  contact  stiffness  k^,  the  plate  flexural 
stiffness  matrix  ,  and  the  impactor  mass  m^.  A  sxmnnary  of  this  study, 

along  with  the  computed  values  of  the  generalized  strains  e  and  e  and  the 

X  y 

contact  parameter  X^,  is  presented  in  Table  I, 
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TABLE  I 


PARAMETRIC  STUDY 

PLATE  IMPACT  CASES  CALCULATED  BY  THE  TIMOSHENKO  SOLUTION  (n=0.459) 


Generalized  Strain 
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i.o«~ 

0.1 


Mass  Ratio,  M  -  inj^/in2 

)C  Experiment  -  Gr/Ep  plate 
— — —  One-degree-of-f reedom  model  ,  Eq«(22) 
_Q - Timoshenko  solution 


Figure  6.  Generalized  strain  curves  for  simply  supported 

rectangular  orthotropic  plates  (rpO.459)  subjected 
to  central  transverse  impact. 
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Note  that  for  all  of  these  calculations  the  value  of  the  aspect-orthotropy 

ratio  ri  is  the  same  (ti  =  0.459).  An  inspection  of  this  table  indicates  that, 

for  any  given  value  of  the  mass  ratio  M,  the  maximum  difference  between  any 

two  values  of  e  or  of  e  is  only  12%,  whereas  the  largest  value  of  X  is  22.6 
X  y  p 

times  the  smallest.  It  may  be  concluded,  therefore,  that  the  dependence  of 
the  generalized  strain  or  on  the  contact  parameter  X^  (or  on  any  other 
parameters)  is  weak,  and  that  for  constant  n  a  single  e  vs.  M  curve  gives 
an  acceptable  representation  of  all  Impact  cases.  Furthermore,  since  the 
powerful  Timoshenko  solution,  which  includes  the  effects  of  contact  behavior, 
indicates  that  such  behavior  does  not  significantly  contribute  to  the  impact 
response,  then  the  validity  of  the  one-degree-of-freedom  model,  which  neglects 
contact  effects,  is  thereby  substantiated  for  the  range  of  low  M. 

The  Timoshenko  solution  may  be  used  to  construct  an  impact  design  curve 
by  plotting  calculated  points  on  an  ^  graph  (Fig,  6).  These  calcu¬ 

lated  points  alone  can  be  approximated  by  a  single  curve  which  may  be  used 
as  a  design  curve.  Note  that  such  a  curve  would  not  be  very  different  from 
the  curve  corresponding  to  the  one-degree-of-freedom  model  discussed  in  a 
previous  section.  Also  included  in  Fig,  6  are  results  of  a  series  of  impact 
experiments  performed  on  the  graphite- epoxy  laminated  plate  described  in 
Appendix  C  [ 8  ] , 

C.  Impact  Experiments 

Several  series  of  impact  experiments  were  performed  on  simply  supported 
and  clamped-edged  plates,  including  strain-measurement  experiments  and  impact- 
to-failure  tests.  The  plate  specimens,  made  of  aluminum  or  graphite/epoxy, 
were  Impacted  by  blunt  steel  projectiles  tising  a  drop-test  apparatus.  In  the 
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strain-measurement  experiments,  the  strain  in  the  plate  was  recorded  using 
strain  gages  and  an  oscilloscope.  In  the  impact-to-failure  tests,  each  plate 
was  impacted  at  gradually  increasing  impact  velocities  until  it  failed.  In 
this  section,  these  experiments  are  described  in  detail,  and  the  results  are 
compared  with  the  theoretically  derived  design  curves. 

1.  Specimens 

The  impact  specimens  included  both  aluminum  and  laminated  graphite/epoxy 
plates.  Three  6061-T6  aluminum  plates  were  used  with  both  simply  supported 
and  clamped  boundary  conditions.  The  dimensions  and  mass  of  these  plates 

are:  (1)  125x125x6.35  mm,  0.282  kg;  (2)  250x125x6.35  mm,  0.561  kg;  and  (3) 
375x125x6.35  mm,  0.843  kg.  These  dimensions  and  masses  include  only  the  area  be¬ 
tween  the  supports  of  the  plates;  an  additional  10-mm  margin  was  left  along  each 
edge  to  support  the  plates* 

Graphite-epoxy  specimens,  fabricated  from  Hercules,  Inc.,  type  AS/35Q1-6 
pre- impregnated  tape  at  various  lay-ups,  were  used  as  simply  supported  and 
clamped  plates.  The  dimensions,  mass,  and  lay-up  of  each  plate  are  listed  in 
Table  II.  These  dimensions  and  masses  again  include  only  the  area  between  the  sup¬ 
ports  of.  the  plates;  a  lO^mm  margin  protruded  beyond  the  supports  around  ^ach  plate. 

Elastic  properties  for  the  composite  plates  were  calculated  from  lamina¬ 
tion  theory,  using  the  following  properties  for  each  layer: 


E 

E 

V 

V 

G 

G 


11 

22 

12 

23 

12 

23 


17.7x10^  psi  =  1.220x10^^  N/m^ 
1.3x10^  psi  =  8.96x10^  N/m^ 

Vi3  =  0.3 

0.2 

=  0.55x10^  psi  =  3.79x10^  N/m^ 
0.54x10^  psi  =  3.72x10^  N/m^ 
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TABLE  III.  SHEAR  AND  FLEXURAL  STIFFNESSES  OF 
GRAPHITE/ EPOXY  PLATE  SPECIMENS 


Flexural 

Stiffnesses' 

(N/m) 


Shear 

Stiffnesses' 

(N/m) 


B-series, 

F-series 

H-series/ 

2470. 

19750. 

9970. 

^2 

710. 

5080. 

4710. 

®22 

963. 

8100. 

13300. 

^6 

-73.3 

-26,7 

-529. 

^6 

-73.3 

-26.7 

-267. 

®66 

748. 

5380. 

4960. 

77700. 

154000. 

' 

147000. 

\5 

0. 

0. 

0. 

*55 

78200. 

155000. 

147000. 

4. 


I 
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TABLE  IV:  DESIGN  CURVE  PARAMETERS  FOR  GRAPHITE/EPOXY  PLATE  SPECIMENS 
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The  computed  values  of  the  flexural  and  shear  stiffnesses  of  each  plate  are 

summarized  in  Table  III.  In  all  further  calculations,  the  D.,  and  D.,  terms 

lo  26 

are  neglected j  that  is,  the  entire  laminate  is  treated  as  specially  orthotropic. 

Calculated  values  of  the  parameter  n  and  the  functions  gj^(n)  and  g2(n) 
necessary  to  construct  the  design  curves  are  listed  in  Table  IV  for  each 
graphite  epoxy  plate. 

2«  Strain-Measurement  Impact  Experiments 

The  plates  were  centrally  impacted  by  blunt  (6.35  ,  12.7,  or  25.4  mm 

contact  radius)  steel  cylinders  of  various  masses  at  several  velocities  using 

a  drop  test  apparatus.  The  strains,  e  and  e  were  measured,  in  directions 

X  y 

parallel  to  the  plate  edges,  directly  opposite  the  impact  point,  using  a 
Micro-Measurements,  Inc.,  type  EA-13-062TT-I20  metal-foil  O^-gO®  rosette 
strain  gage,  and  recorded  on  a  Tektronix  model  565  oscilloscope.  Impact  ^ 

velocities  were  calculated  from  the  drop  height.  ' 

For  Impacts  on  clamped  plates,  the  plates  were  clamped  along  all 
edges  by  a  frame  of  rectangular  steel  bars  (Fig.  7).  The  entire  assembly  was 
fastened  to  a  large  lead  plate  and  rested  on  a  concrete  floor. 

In  the  experiments  performed  on  simply  supported  aluminum  plates,  each  plate  i 

was  rested  on  top  of  a  frame  of  rectangular-cross-section  steel  bars.  For  ] 

the  composite  plates,  however,  the  clamping  device  mentioned  above  was  modi-  i 

fied  to  simulate  a  simply  supported  boundary.  Grooves  to  accommodate  3.175  mm  j 

(1/8")  rods  were  machined  in  each  member  of  the  clamping  device  such  that, 
when  the  clamping  device  was  assembled,  the  specimen  plate  would  be  resting 
on  the  rods. 

Experimental  results  of  the  strain-measurement  tests  on  both  clamped 
and  simply  supported  plates  are  presented  in  Appendix  E  in  both  tabular  and 
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Figure  7.  Exploded  viev  of  plate  clamping  device. 
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graphical  (e  vs.  M)  form.  Also  Included  in  Appendix  E  are  some  typical  strain 
vs.. time  oscilloscope  traces.  Figures  8  and  9  are  typical  examples  of  the  e  vs. 

M  curves  for  the  clamped  plates  and  the  simply  supported  plates,  respectively. 

In  each  of  these  figures  a  dashed  line  has  been  fitted  by  eye  through  the 
experimental  points.  The  maximum  percentage  difference  (variation)  between 
this  line  and  the  experimental  data  is  also  shown  on  each  figure.  For  both 
the  clamped  and  the  simply  supported  cases,  the  largest  value  of  this  variation 
is  36%. 

Experimental  data  from  specimens  having  approximately  the  same  value  of 
the  parameter  n  have  been  presented  on  the  same  design  curve  except  the  three 
smallest  sized  plates _ (Bl, FI,  and  HI).  Values  of  e  for  these  three  plates 
are  considerably  smaller  than  those  for  larger  plates  with  the  same  parameter 
ri.  It  is  believed  that,  due  to  the  small  size  of  these  plates,  the  boundary 
constraint  has  more  influence,  and  the  deflection  mode  is  different  from  the  other 
plates.  Perhaps  thick-plate,  damping,  and  shear  effects  are  more  pronounced. 

Data  from  these  plates,  although  presented  on  the  e  vs .  M  plot  of  Figs.  10  and 
11,  were  not  used  in  constructing  the  experimental  design  curve  line  nor  in 
calculating  the  variation. 

Also  shown  in  Figs.  8  and  9,  and  in  each  design  curve  in  Appendix  E,  is 
the  theoretical  (one-degree— of -freedom)  design  curve.  Comparison  of  this 
curve  and  the  one  from  experimental  data  shows  an  average  discrepancy  of  about 
40%.  Note  that  the  theoretical  curve  is  always  conservative  (l.e.,  predicts 
a  higher  plate  strain)  over  the  range  of  experimental  results.  However,  the 
slope  of  the  line  drawn  through  the  experimental  data  is  in  all  cases  roughly 
the  same  as  the  slope  of  the  theoretical  curve. 
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Plgare  8,  Typical  M  Curve  for  Clamped  Orthotropic  Elates 
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Figure  11.  Strain  Curve  for  Simply  Supported  Plates 
Impacted  by  Large  Masses 
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previously  discussed,  ell  inipect  dete  cen  be  plotted  on  s  singls 
curve,  T  vs.  M.  This  is  done  for  impact  data  on  claii?)ed  plates  in  Fig.  10 
and  on  simply  supported  plates  in  Fig.  11.  Note  that,  since  £  is  independent 
of  both  the  parameter  ri  and  the  boundary  conditions,  both  cases  could  have  been 
plotted  on  a  single  figure.  Because  of  the  overlapping  of  experimental  data, 
the  clamped  and  simply  supported  cases,  were  presented  separately  for  clarity. 

3.  Impact- to~Failure  Experiments 

In  addition  to  the  strain  measurement  impact  experiments,  a  series  of 
tests  were  conducted  in  which  laminated  composite  platds  were  impacted  to 
failure.  In  each  of  these  tests,  a  plate  was  repeatedly  impacted  by  the 
same  projectile  using  a  drop— weight  apparatus  at  gradually  increasing  veloci¬ 
ties  (drop  heights)  until  failure  of  the  specimen  was  detected.  The  specimen 
was  examined  for  failure  after  each  impact  both  visually  and  ultrasonically, 
using  a  hand  held  pulse-echo  transducer  and  reflectoscope. 

The  average  between  the  highest  impact  velocity,  v^»  for  which  a  particular 
specimen  did  not  fall,  and  the  lowest  velocity,  v^,  for  which  any  failure  was 
detected,  is  regarded  as  the  critical  or  failure  velocity,  v^,  or 


This  velocity  is  used  in  characterizing  the  impact  resistance  of  each  specimen. 

The  results  of  these  experiments  are  presented  in  both  graphical  and 
tabular  form  in  Appendix  E.  Figure  12  is  a  typical  example  of  the  v^vs.  M 
curves  for  the  clamped  plate.  Again,  a  dashed  line  has  been  drawn  to  indicate 
the  design  curve  based  only  on  experimental  results.  As  can  be  seen,  there 
is  good  agreement  between  the  experimental  data  and  the  theoretical  one- 
degree-of-freedom  curve. 

Again,  a  design  curve  Independent  of  the  parameter  n  and  the  boundary 
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conditions  (v  vs.  M)  may  be  used  to  present  Impact  data.  However,  from 
Eq.  (34).  it  is  apparent  that  only  those  plates  with  the  same  value  of 
should  be  presented  on  this  curve.  For  the  composite  plates,  therefore,  a 
V  vs.  M  curve  should  be  prepared  for  each  different  lay-up.  Figure  13  is  i 
such  a  plot  for  the  experiments  performed  on  B-serles  plates. 

After  damage  was  detected,  several  plates  were  ultras onically  C-scanned. 
Transducer  output  was  passed  through  an  analog  to  digital  converter.  Data 
was  then  processed  by  Fourier  transform  techniques  using  the  first  five  wave¬ 
forms.  Typical  transducer  response  is  shown  in  Fig.  14  for  an  undamaged 
area  of  a  plate.  The  top-surface  and  bottom-surface  echoes  are  clearly  visible 
with  little  evidence  of  internal  reflection.  Figures  15  and  16  show  trans¬ 
ducer  response  over  damaged  areas  of  two  different  specimens.  Figure  15 
shows  a  damaged  region  extending  from  about  one-quarter  of  the  plate  thick¬ 
ness  below  the  surface  to  the  midplane  of  the  plate.  This  type  of  damage 
has  been  observed  by  use  of  the  hand-held  transducer  in  approximately  75%  of 
the  failure  tests.  It  is  believed  that  this  mode  of  failure  is  delamlnatlon. 
Figure  16  is  representative  of  the  remaining  25%  of  the  failure  tests.  As 
can  be  seen,  there  is  little  Indication  of  damage  through  the  thickness  of 
the  plate.  This,  along  with  the  absence  of  a  bottom-surface  echo  indicates 
that  the  plate  has  probably  been  damaged  near,  the  top  surface.  It  is  believed 

that  this  mode  of  failure  results  from  fiber  breakage,  due  to  bending  and/or 
contact  effects. 

-Design  Curves  for  Shear  Failure  Due  to  Impar^ 

A  structure  subjected  to  impact  may  also  fall  due  to  the  high  level  of 
shear  stress  produced.  In  this  section,  a  theory  based  on  the  one-degree-of- 
freedom  Impact  model  is  developed  for  predicting  the  maximum  impact-induced 

transverse  shear  stress  in  simply  supported  beams  and  plates. 

40 


NADC-78259-60 


Typical  Oritloal  Tulooity  Ounre,  ^vs.  M,  for  CXaijed 
Orthotropic  Plates:  Esiperimental  Data  from  Specimens 
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Figure  13.  Critical  Impact  Velocity  Curve  for 
Clamped  and  Simply  Supported  Plates 
Impacted  by  Large  Masses. 
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Figure  16.  Typical  Reflectoscope  trace  for  plate  specimen  with  top-surface  damage. 
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1.  Simply  Supported  Beam 

By  extending  the  one-degree-of-freedom,  energy-conserved  model  previously 
discussed,  a  design  curve  for  predicting  the  peak  shear  stress  in  a  simply 
supported  beam  subjected  to  transverse  liiq>act  may  be  generated. 

the  maximum  shear  stress  in  the  beam  is  related  to  the  maximum 
deflection  by  the  static  beam  equation.  The  deflection  of  a  simply  supported 
beam  subjected  to  a  static  load  at  midspan  Is 


w(x)  =  w 


1  L 


0<x<L/2 


(54) 


where  w^  is  the  deflection  at  midspan.  The  shear  force  generated  is 


V  =  El 


d\ 

dx3 


O/  El 
24^w^ 


(55) 


For  a  beam  of  rectangular  cross-section,  the  maximum  shear  stress  occurs  at 
the  neutral  axis  and  is  given  by 


/ 


T 

max  2A 

36EI  (56) 

- 3  w 

AL-^  mx 

The  maximum  deflection  may  be  related  to  the  impact  velocity  by 

■‘‘max 

the  conservation-of-energy  condition,  which  implies  that  the  initial  kinetic 
energy  of  the  Impactor  is  entirely  converted  into  bending— strain  energy  in 
the  beam. 


46 


NADC-78259^ 


By  defining  the  dimensionless  (generalized)  shear  stress  as 


T  = 


T  2. 

T  L  a  A 
max _ 

vEI 


(58) 


and  combining  Eqs.  (56)  and  (57),  we  obtain  the  eqijation  for  the  shaar-stress 


design  curve. 

'-31^  (59) 

This  eqiiatlon  Is  plotted  in  t  vs.  M  coordinates  In  Fig.  17. 


a.  Relative  Importance  of  Shear  and  Bending 

Depending  on  geometry,  strength,  and  elastic  properties,  a  particular 
beam  subjected  to  low-velocity  Impact  may  experience  failure  Initiated  by 
bending  or  shear  effects.  The  theories  presented  here  for  predicting  bending 
strain  and  shear  stress  due  to  Impact  are  both  based  on  the  static  beam 
deflection  curve.  Therefore,  the  relative  Importance  of  bending  and  shear 
effects  may  be  estimated  by  considering  the  static  relations  between  shear 
and  bending. 

If  a  maximum-stress  failure  criterion  Is  adopted  for  both  shear  and 
bending,  then  failure  due  to  shear  will  be  more  likely  If,  for  any  Impact 
velocity  V,  we  have 


-2- 


(60) 


‘f  “f 

where  the  subscript  f  denotes  values  at  failure.  By  combining  this  relation 


with  Eqs.  (12),  (13,  (59),  and  (50),  we  obtain 

2 

LXf 


(61) 


Based  on  the  approximate  theory  presented  here,  this  Inequality  must  be 
satisfied  for  Impact  failure  to  Initiate  by  shear  effects;  If  It  Is  not 
satisfied,  failure  will  Initiate  by  bending. 
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b.  Critical  Impact  Velocity  Curve 

The  relative  Importance  of  shear  and  bending  effects  In  Initiating 
failure  of  a  particular  beam  may  be  Illustrated  by  presenting  the  design 
curve  In  a  form  similar  to  that  presented  In  Section  II  for  plate  Impacts. 

In  this  alternate  form,  theoretical  design  curves  and  experimental  data  are 
plotted  In  coordinates  of  a  dimensionless  Impact  velocity  v  or  v  and  the  mass 
ratio  M.  Then,  each  critical  velocity  curve  divides  the  plane  Into  two  regions; 
above  the  critical  bending  curve  Is  the  region  bending  failure  occurs,  and 
below  Is  the  no-bendlng-fallure  region.  The  same  Is  true  for  the  critical 
shear  curve. 


Specifically,  If  the  dimensionless  Impact  velocity  for  the  beam  Is  defined 


as 


V  = 


vh 


(62) 


a 


then  the  one-degree-of-freedom  energy-conserved  Impact  model,  predicts  that 
failure  will  Initiate  In  the  beam  due  to  bending  for 


2e 

V  =  -7=^ 

c  ^ 


(63) 


where  Is  the  tensile  failure  strain,  assuming  a  maximum— strain  failure 
criterion. 

Similarly,  according  to  the  theory  discussed  above  for  estimating  the 
shear  stress  due  to  Impact,  Eq.  (59),  failure  due  to  shear  occurs  for 


V  =  -= -  ✓S 

c  /3  Eh 


(64) 


where  Is  the  shear  stress  at  failure,  assuming  a  maximum— stress  failure 


criterion  In  shear. 
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Equations  (63)  and  (64),  plotted  in  v  vs.  M  coordinates,  represent  the  two 
critical  velocity  design  curves.  The  relative  importance  of  shear  and  bending 
regarding  failure  of  a  particular  structure  is  immediately  apparent  from  such 
a  plot  —  the  lower  curve  predicts  failure  at  a  smaller  impact  velocity,  indi¬ 
cating  that  the  failure  mode  (shear  or  bending)  associated  with  the  lower 
curve  is  the  critical  one. 

Consider,  for  example,  the  impact  failure  tests  on  22  graphite-epoxy  beams. 
These  beams  were  fabricated  from  Hercules  AS3501  at  a  lay-up  of  [+  45/O2/+  45]2» 
All  beams  have  a  span  to  thickness  ratio  (L/h)  of  48.53.  Data  from  all  the 
impact  tests  (with  or  without  failure)  performed  on  these  beam  specimens  are 
tabulated  in  Appendix  E  and  plotted  in  v  vs.  M  coordinates  in  Fig.  18.  Also 
included  in  this  figure  are  curves  representing  Eqs.  (63)  and  (64).  The  value 
of  the  bending  failure  strain  used  is  0.0168;  the  value  of  the  shear  failure 
stress  used  is  the  value  of  the  interlaminar  shear  strength  given  by  the 
material  manufacturer  (Hercules  Product  Data  Sheet  No.  832)  as  18,900  psi 
(130.  MN/m^). 

Note  that  the  lower  curve  ^n  Fig.  18,  which  represents  Eq.  (63),  approxi¬ 
mately  divides  the  experimental  points  according  to  whether  the  points  correspond 
to  failure  or  no-failure  tests.  Also,  the  curve  corresponding  to  shear  failure 
falls  much  higher  on  this  plot  than  the  curve  for  bending  failure;  this  indicates 
that  for  these  beams  the  bending  is  the  mode  in  which  failure  initiates. 

2.  Simply  Supported  Plate 

A  design  curve  for  estimating  the  peak  transverse  shear  stress  at  the 
edge  of  a  siii5>ly  supported,  rectangular,  orthotropic  plate  subjected  to  central 
lateral  impact  may  also  be  developed. 

It  must  be  recognized  that  the  edge  shear  is  only  a  rough  estimate  of 
the  maximum  shear  stress  occurring  in  the  plate;  much  higher  stresses  may  be 

generated  in  the  immediate  vicinity  of  the  Impact  point.  However,  analysis 
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of  the  problem  by  plate  theory  alone  (as  is  done  here)  leads  to  a  singularity 
in  the  shear  at  this  point,  due  to  treatment  of  the  impact  load  as  a  concen¬ 
trated  force.  A  more  sophisticated  analysis  would  focus  on  the  three-dimen¬ 
sional  stress  field  near  the  impact  location.  Still,  for  design  purposes, 
knowledge  of  the  edge  shear  can  sometimes  be  useful.  If,  for  a  particular 
impact  situation,  it  exceeds  the  allowable  shear  stress,  failure  will  certainly 
occur;  on  the  other  hand,  a  low  edge  shear  stress  does  not,  of  course, 
guarantee  survival  of  the  structure. 

As  in  the  previous  case  of  the  beam,  the  maximum*  edge  shear  stress  in 
the  plate  is  first  related  to  the  maximum  central  deflection  by  assuming  that 
the  static  relations  hold.  In  terms  of  generalized  plate  forces,  the  maximum 
values  of  the  transverse  shear  stress  components  are 


xz  2h 

V.  ■  2h  (^5) 

for  a  homogeneous  plate.  For  a  specially  ortho tropic  plate,  the  generalized 
forces  are  related  to  the  deflection  distribution  by 


Q  =  -  /  D 

X  9x  't 


-L/d  + 

V  11  9,2  3^2 


where 


0  =  -  H  —  +  D  — ^ 


“  =  “12  +  2D55 


It  is  again  assumed  that  the  dynamic  deflection  due  to  impact  is  the 
same  as  the  static  deflection  due  to  a  central  point  load.  In  this  case,  the 
maximum  edge  shear  force  occurs  at  the  middle  of  the  edge.  If  it  is  also  assumed 


that  H  =  "IP 2^2^22’ 
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Figure  19.  Edge  shear  stress  curves  for  simply  supported  rectangular  orthotropic  plates. 
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III.  IMPACTS  BY  SMALL  IMPACTORS 


In  this  section,  approximate  equations  for  predicting  the  peak  strain 
response  of  simply  supported  beams  and  plates  to  impacts  by  small  impactors 
are  derived.  The  accuracy  of  these  equations  is  demonstrated  by  comparison 
with  experimental  results  and  with  the  Timoshenko  solution. 

A.  Simply  Supported  Beam 
1.  Approximate  Analysis 

It  has  been  observed  in  experiments  and  in  calculations  of  the 
Timoshenko  solution  that  for  the  case  of  small  impactors  (mass  ratio  M 
greater  than  about  2),  the  interaction  between  beam  and  impactor  is  a  single, 
sudden  blow  of  short  duration  compared  to,  say,  the  fundamental  period  of 
vibration  of  the  beam.  This  is  quite  different  from  the  large  impactor  case 
where  multiple  impacts  occur.  By  assuming  that  the  beam  does  not  appreciably 
deflect  during  the  short  period  of  contact  and  by  making  a  few  other  approxi¬ 
mations,  we  may  derive  a  relation  which  is  useful  in  generating  design  curves. 

The  assumption  that  the  beam  does  not  deflect  during  contact  is  equiva¬ 
lent  to  assuming  that  the  impactor  is  much  less  massive  than  the  beam,  which 
may  thus  be  treated  as  a  semi-infinite  body.  An  approximate  solution  for  the 
contact  force  during  the  elastic  impact  of  a  spherical  impactor  against  a 
flat  semi- infinite  body  has  been  calculated  by  Hunter  [9]  as 


where 


_  |p  sin  5t  0  <  t  <  Tf/g 
1^0  t  >  tt/C 


P  =  1.068  V  m2  S 


5  =  0.9768 


2.1/5 
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The  response  of  a  simply  supported  beam  to  a  centrally  applied  trans¬ 
verse  dynamic  load  of  the  general  form  of  Eq.  (74)  is  derived  in  Appendix  D 
The  maximum  deflection  and  bending  strain  occur  at  midspan  at  a  time  about 
one-quarter  of  the  fundamental  period  after  impact.  The  peak  strain  may  be 
expressed  as 


e 


max 


PLh 

El 


f(6) 


(77) 


where 


and 


f(6)  i 

IT  i=l,3,5 


COS 


1  fundamental  period  of  beam 

2  contact  force  duration 


) 


(78) 


(79) 


Substitution  of  the  approximate  value  of  5  given  in  Eq.  (76)  into 


Eq.  (79)  yields  an  approximate  value  for  g. 

2,1/5 


g  =  (0. 09897)' 


/L 


2  V  /vk, 


(80) 


This  quantity  may  be  expressed  in  terms  of  the  dimensionless  parameters 
previously  derived  as 


2x1/5 


g  =  (0.1564) (  y-  j 


(81) 


b 


where  M  and  are  defined  in  Eqs.  (3)  and  (53),  respectively. 
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By  substituting  for  P  and  3  into  Eq.  (77)  according  to  Eqs.  (75)  and 
(76),  and  recalling  the  definition  of  generalized  strain,  Eq.  (12) 
we  obtain 

e  =  1.068  TT^  f  (  3  )  ^  (82) 

This  result  may  Itself  be  used  to  generate  a  series  of  design  curves  by 
plotting  e  vs.  M  for  different  values  of  3  ,  however,  this  equation  may  be 
further  simplified  as  demonstrated  below. 

Note  that  Eq.  (82)  shows  a  strong  dependence  of  the  generalized  strain 
e  on  both  the  mass  ratio  M  and  the  parameter  3  .  Also,  observe  that  for  a 
constant  value  of  3  ,  e  is  inversely  proportional  to  M,  rather  than 
as  for  the  case  of  large  impactors.  This  suggests  that  for  small  impactors 
the  peak  strain  is  related  to  the  momentum  of  the  impactor,  rather  than  to 
its  kinetic  energy  as  for  large  impactors.  That  is,  large-M  cases  are  governed 
by  the  impulse  exerted  by  the  impactor,  whereas  small-M  cases  are  character- 
ized  by  the  work  done  on  the  beam. 

The  influence  of  the  value  of  the  parameter  3  on  the  generalized 
strain  s  can  be  estimated  by  taking  the  limit  of  Eq.  (82)  as  3  becomes  large. 

In  Table  V,  calculated  values  of  f(3  )are  presented  (see  Appendix  D). 
Examination  of  this  table  reveals  that,  as  3  increases,  the  quantity  /3  f(3) 
approaches  the  value  0.12412.  Observe  also  that,  for  3  >  7,  /3  f(3)  never 
differs  from  this  value  by  more  than  2%.  Since  most  practical  problems  satisfy 
the  condition  3  >  7,  we  can  assume  with  little  inaccuracy  that 


/F  f(  e )  ^  0.1241  . 


NADC -78259 .60 


TABLE  V 

COrn'IER-CALCULATED  VALUES  OF  f(8) 


f(B) 

^  f(6) 

3.0 

0.06466 

0.11200 

7.0 

0.04590 

0.12145 

10.0 

0.039C3 

0.12342 

15.0 

0.03190 

0.12356 

20.0 

0.02780 

0.12432 

30,0 

0.02264 

0.12402 

50.0 

0.01755 

0.12408 

80.0 

0.01388 

0.12414 

3  00,0 

0.012413 

0.12413 

150.0 

0.010132 

0.12412 

200.0 

0.008777 

0.12413 

300.0 

0,007166 

0.12412 

500.0 

O.CO:)551 

0.12412 
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Substitution  into  Eq.  (82)  yields 


e  ^  1.308 

/r“  M 


(83) 


or,  defining  a  new  generalized  strain  for  the  small  impactor  case  as 


e*  E 


/T 


(84) 


we  then  have 


E*  = 


1.308 

M 


(85) 


This  simplified  equation  may  be  used  as  a  design  curve  for  the  case  of  small 


impactors  by  plotting  e*  vs.  M.  In  Fig.  21,  such  a  curve  is  compared  with 
the  experimental  data  recorded  in  Appendix  E.  In  general  the  experimentally 
determined  values  of  e*  are  roughly  60%  lower  than  those  predicted  by  Eq.  (85). 
Note,  however,  that  the  band  of  experimental  points  on  this  graph  form  a 
relatively  narrow  band;  for  any  value  of  M,  the  variation  in  e*  is  at  most 
i  25/  (et  about  M  —  6)  and  much  less  for  the  higher  range  of  M  (greater  than 
about  10).  This  variation  in  e*  is  much  smaller  than  the  variation  in  e 
^lone,  for  the  same  data.  Also,  the  slope  of  this  band  is  very  close  to 
that  of  the  curve  corresponding  to  Eq. (85) .  This  suggests  that  the  locus 
of  the  experimental  data  plotted  in  these  particular  coordinates  may  itself 
be  useful  as  a  design  curve.  We  may  conclude  that  Eq.  (85),  while  not  accurate, 
IS  significant  in  that  it  serves  to  identify  the  form  of  the  relationship  among 
the  several  parameters  governing  beam  response  to  impacts  by  small  masses. 

2 .  Timoshenko  Solution 

Also  Included  in  Fig.  21  are  several  points  corresponding  to  calcula- 
bions  using  the  Timoshenko  solution;  data  from  these  calculations  are  summar¬ 
ized  in  Table  VI.  Note  that  these  points  form  a  very  narrow  band,  the 
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TARLE  VI 


CALCULATIONS  USING  TIMOSHENKO 


SOLUTION  FOR  LARGE  ll\SS  1L\TI0  M 


Beam  Mass, 
Dimensions, 
and  Material 


0.389  kg 

197  X  16  X  16  Iran 


steel 


Mass  Impact  Co 

Ratio,  Velocity,  Sti 

M=m, /m„  V  (m/sec)  k«( 


S.72  2.466x10 


Generalized 

Strain, 

—  2 , 

£  =  ea  /hv 


0.1807 


0.1894 


0.2835 


0.5727 


0.5865 


2.415x10^^'  0,1782  '  0.0416 


0.0405 


0,0412 


0.0734 


0.0728 


0.0757 


0.0720 


0.0744 


0.1287 


0.2013 


0.1984 


0.389  kg,  steel 
394  X  8  X  16  mm 


0.389  kg,  steel 
394  X  16  X  16  mm 


0.389  kg,  steel 
197  X  8  X  32  mm 


0.1943  kg,  steel 
197  X  8  X  16  r.m 


13.84 


8.72  2.466x10 


0.2502 


0.0322 


.2703  0,0361 


0.3172  0.0559 


0.0789 


0.1192 
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variation  being  less  than  +  5%;  again,  the  variation  in  e*  for  the  same 
data  is  much  larger*  In  addition,  the  band  of  points  is  much  closer  to  the 
curve  corresponding  to  Eq.  (85)  than  are  the  experimental  points;  these  points 
lie  only  15  to  45%  below  the  curve  and  appear  to  approach  the  curve  as  M 
increases.  The  closeness  of  these  results  Implies  that  the  lack  of  agree¬ 
ment  between  Eq.  (85)  and  the  experimental  data  is  not  only  due  to  the  approx¬ 
imations  used  in  deriving  the  equation.  Apparently,  there  are  additional 
effects  which  significantly  influence  beam  response  to  impacts  by  small  masses 
which  are  hot  accounted  for  even  in  the  powerful  Timoshenko  solution.  Several 
researchers  have  modified  the  basic  solution  to  include  some  of  these  effects, 
but  these  more  involved  models  are  too  cumbersome  for  the  development  of  an 
impact  design  curve,  since  for  each  effect  considered,  another  parameter  (or 
possibly  several)  must  be  introduced.  Such  complications  are  surely  justified 
when  one  is  analyzing  a  specific  Impact  problem,  but  are  of  little  assistance 
when  one  is  designing  a  structure  to  resist  potential  Impacts  due  to  a  range 
of  masses,  velocities,  materials,  etc. 

B.  Simply  Supported  Plate 

In  this  section  is  developed  a  design  curve  for  predicting  the  response 
of  a  simply  supported  plate  to  impact  by  a  small  mass  using  the  same  approach 
as  applied  above  to  a  beam. 

1.  Approximate  Analysis 

As  in  the  case  of  the  beam,  we  assume  that  the  plate  does  not  appreciably 
deflect  during  the  short  period  when  the  Impactor  is  actually  in  contact  with 
the  plate.  Again,  this  assumption  is  equivalent  to  supposing  that  the  impactor 
is  much  less  massive  than  the  plate.  Thus,  the  contact  force  is  roughly  given 
by  the  approximate  solution  for  a  sphere  elastically  striking  a  flat  semi¬ 


infinite  body,  Eqs.  (74),  (75),  and  (76). 
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The  dynamic  response  of  a  simply  supported  plate  to  a  central  transverse 
load  F(t)  is  given  by  [  6 ] 


,  f  .^(m+n-2)/2 

W  =  —  I  F  -  sin  —  sin  SIE2. 

m,  f-  ^  0)  a  h 


-t  60 

1  m  n  nm 

t 


£ 


F(t)  sin  (1)  (t-T)dT 
mn 


(86) 


For  a  thin  specially  orthotropic  plate,  the  natural  frequencies  are  given 

by  [10] 

4 


2  TT 
60  =  — 
mn  ph 

If  we  again  assume  that 

then  we  have 

in  which 


4  2  2  4 

D  £ — h  2fD  +2D  )  £■  Jl +  T) 

^11  4  ^  2^2  ^  ®22  ,4 

a  abb 


"^12  (43) 


2  ^  2 
0,  =  co„  2-+J3n_ 

mn  11  1+n 


(87) 


1/2 


«+’>) 

3 


(88) 


Substitution  of  the  assumed  form  of  the  contact  force,  Eq.  (74),  into 
Eq.  (86)  and  integration  yields 

m+n-2 

^  sin  SS  sin  ^ 


w  - 


1 1 


l+n 


m,  0), ,  ^  2  2 

1  11  m  n  m  +  nn 


2P 


(-1) 


5(l-m  /C) 

mn 


y — 5—  sin  0)  (t+ir/O  cos 
mn 


0)  tr 
mn 


2? 


(89) 


In  the  corresponding  formula  for  the  beam,  all  terms  in  the  series  reach 
their  maximum  values  at  exactly  the  same  time,  due  to  a  simple  relationship 
which  exists  among  the  natural  frequencies  (see  Appendix  D,  Eq.  (Dl)). 
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Now,  if  we  substitute  for  P  in  Eq.  (92)  according  to  Eq.  (75),  we 
obtain 


e  =  4(1.068)hv 
x,inax  |°11 

3  ir  A 

p  mn 

mn  ^ 

Finally,  If  we  define  the  generalized  strain  for  impact  of  a 
plate  by  a  small  mass  as 


n  - 

m  n  mn  3 


(97) 

simply  supported 


*  _  x.max _ _ 

X  4(1.068)hv  3  g(3  ,ri) 


(98) 


where 


B 


7T  A 


“  l+n  ^  ^  f~  — 2~^ - 2  2B"”” 

P  m  n  3  -  A  ^  ^3. 


(99) 


p  mn 

then  the  equation  for  the  design  curve  according  to  this  approximate  solution 
is 

e*  =  l/M  (100) 


Note  that  the  function  g(3p,q),  plotted  in  Fig.  22,  is  practically 

independent  of  the  value  of  3  ,  at  least  for  realistic  values  of  3  (3  >7). 

P  p  p 

In  Fig.  22a,  g(3p,n)  is  plotted  over  n  for  two  widely  different  values  of 

3p,7  and  500;  the  two  curves  are  almost  the  same.  In  Fig.  22b,  the  maximum 

and  minimum  values  of  g(3„,n)  for  3  in  the  range  7  <  3  <  500  are  plotted 

P  P  p  ““ 

for  each  value  of  rj;  again,  the  two  curves  are  quite  close,  indicating  that 

g(6„»Tl)  does  not  vary  much  with  3  . 

P  p 
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Figure  23.  Strain  Curve  for  Impact  of  Simply  Supported 
Rectangular  Plate  by  a  Small  Mass. 
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2.  Comparison  with  Experimental  Results 

In  Fig.  23,  the  design  curve  for  impact  of  plates  by  small  masses  is 
presented  along  with  the  results  from  experiments  performed  at  Drexel.  Two 
graphite-epoxy  plates  which  were  also  used  in  the  large  impactor  experiments, 

FI  and  F2,  were  impacted  by  aluminum  projectiles.  These  experiments  were 
performed  using  an  air  gun  with  the  projectile  velocity  recorded  by  a  photo¬ 
diode  system.  Data  from  these  experiments  are  tabulated  in  Appendix  E.  As 
can  be  seen  in  Fig.  23,  the  agreement  between  the  theoretical  design  curve  and 
the  design  curve  based  on  experimental  results  is  of  the  same  order  as  in  the 
large  impactor  case;  that  is,  the  experimental  cunre  is  about  40%  below  the 
theoretical  curve  and  has  approximately  the  same  slope. 

Also  tabulated  in  Appendix  E  are  experimental  data  found  in  Schwieger 
111].  These  experiments  were  performed  on  a  square  Duraluminum  (AlC’iMgl) 
plate,  550  x  550  x  4.97  mm,  centrally  impacted  by  steel  spheres  of  20--  and 
30-  mm  radius.  Tliese  data,  when  plotted  in  e*  vs.  M  coordinates,  do  not  corapar 
well  with  either  the  Drexel  experiments  or  the' ’theory.  One  explanation  of  this 
discrepancy  lies  in  the  size  of  the  plate  tested  by  Schwieger.  In  the  Drexel 
experiments,  the  ratio  of  the  plate  span  to  thickness  (L/h)  has  a  value  of  18 
whereas  in  the  e:q?eriments  performed  by  Schwieger,  this  ratio  has  a  value  of  11 
Thus,  according  to  Zener  [12],  the  response  of  the  plate  used  by  Schwieger  is 
governed  by  membrane  and  wave  propagation  effects.  Since  these  effects  are 
not  treated  by  the  current  theory,  the  data  from  Schw^ieger  are  not  included 


on  the  design  curve. 
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IV.  DESIGN  PROCEDURE 


In  this  section,  procedures  are  recommended  for  using  the  design  curves 
presented  in  this  report  and  for  constructing  new  design  curves  for  structures 
not  treated  here  based  either  on  analytical  calculations  or  on  experimental 
data.  It  is  emphasized  that  these  design  curves  are  limited  to  predicting 
only  the  maximum  structural  response  to  low-velocity  impact. 

A.  Use  of  Design  Curves  in  This  Report 

The  design  curves  in  this  report  predict  the  peak  bending  strain  in  simply 
supported  and  clamped  plates  impacted  by  large  impactors,  and  in  simply  supported 
beams  and  plates  impacted  by  small  impactors.  Also,  in  [1],  simply  supported 
beams  impacted  by  large  masses  were  treated  in  detail.  Methods  for  the  use  of 
each  of  these  design  curves  will  be  discussed  separately  here. 

1.  Simply  Supported  Beam  ^ 

Design  curves  for  simply  supported  beams  impacted  by  large  Impactors  and 
by  small  impactors  are  repeated  here  for  convenience  (Figs.  24  and  25,  respec¬ 
tively).  In  each  case,  we  have  selected  a  design  curve  fit  by  eye  to  the 
experimental  results  as  the  one  we  will  use. 

The  first  step  in  using  these  curves  is  to  compute  the  mass  ratio 

„  ™1  structural  mass  _ _ _ 

M  =  —  =  — : - - -  (101) 

m2  impactor  mass  '  ' 


1 


If  M  is  less  than  2,  the  impact  case  falls  into  the  domain  of  large 
impactors;  if  M  >  2,  small  impactors. 
a.  Large  Impactors 

From  Fig.  24,  read  the  value  of  the  generalized  strain  e  corre¬ 
sponding  to  the  value  of  M.  Then,  using  this  value  e,  compute  the  tnflv^1ll11Tn 
Impact  strain, 

I 

i 

\ 
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_ SELECTED  DESIGN 

‘curve 


MASS  RATIO,  M=mj/in2 


Figure  25. 


Selected  design  curve  for  Impacts  of  simply  supported 
beams  by  small  Impactors. 
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e  =1^  (102) 

max  kc 

wliGTT©  ^  ~  Idgsih  d.6ptli 

k  =  radius  of  gyration  of  beam 
cross-section 

V  =  impact  velocity 

c  =  i/e/p  =  speed  of  flexural  waves 
in  the  beam 

The  use  of  this  curve  will  be  illustrated  by  an  example. 


Example  #1 

3 

For  a  simply  supported  steel  (density,  7.9  g/cm  )  beam  of  rectangular 
cross-section  with  dimensions  L  x  h  x  b  =  500x25x20  mm  inq)acted  at  10  m/s 


by  a  mass  of  4.0  kg,  we  have 

„  “l  L  bh  _  (7900) (.50) (0.025) (0.020) 

4.0 

-  0.494 

From  Fig.  24,  the  value  of  e  corresponding  to  M  =  0.494  is  e  =  1.20.  Then, 
by  Eq.  (102),  we  have 


max 


=  (1.20) 


(0.025) (10) 


0.02^  210x10^ 


4/2 


\ 


7900 


=  0.0081 


b.  Small  Impactors 

From  Fig.  25,  read  the  value  of  the  generalized  strain  e*  corres¬ 
ponding  to  the  value  of  M.  Then  compute  the  maximum  of  impact  strain. 


-  /g  is 


where 


and 


3  =  (0.09897) 


ck 


(103) 


-  Hertzian  contact  stiffness. 
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Example  //2 

Consider  the  beam  of  Example  #1  impacted  by  a  steel  sphere  of  radius 
20  mm  at  20  m/s.  Then,  by  Eq.  (101),  we  have 

M  =  =  (7900) (0.50) (0.025) (0.020) 

“2  ir  (0.020)^  (7900) 

=  7.46 

From  Fig.  25,  the  value  of  e*  corresponding  to  M  =  7.46  is  e*  =  0.060.  The 
Hertzian  contact  stiffness  is 


2  Stt 


Et7  A  V  Eir  J. 

beam  \  ''impactor_ 


=  /0.020 

jir 


[— 

l-ir  21 


210x10' 


so  that 


=  2.176x10^°  N/m^^^ 


(0.50] 


3  =  (0.09897) 


=  28.1 

Therefore,  the  maximum  bending  strain  is 


1/210x10^ 

0.025 

f  7900 

strain  is 

2 

20(2.176x10^^) 

0.265^ 


e  =  (0.060)  /2^ 
max 


(0.025) (20) 

9 


210x10^  (0.025) 
'  7900  /J2 


=  .00427 
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2.  Simply  Supported  and  Clamped  Plates 

The  procedures  for  using  the  design  curves  for  rectangular  plates  is 
similar  to  those  outlined  above  for  the  beam.  For  convenience,  we  have  re¬ 
peated  here  the  design  curves  for  simply  supported  and  clamped  plates  impacted 
by  large  impactors  (Figs.  26  and  27),  and  for  simply  supported  plates  impacted 
by  small  Impactors  (Fig.  28).  The  first  step  in  using  any  of  these  curves  is 
to  compute  the  mass  ratio  M  according  to  Eq.  (101)  and  the  aspect— orthotropy 


ratio. 


”  (D  #22^^11 


(104) 


a  and  b  are  the  plate  dimensions  in  the  x  and  y  directions.  For  the 
plates  treated  here,  impact  cases  in  which  M  <  4  fall  into  the  domain  of 
large  impactors,  while  cases  where  M  >  4  are  in  the  domain  of  small  impactors. 
The  use  of  the  design  curves  in  each  of  these  domains  will  be  discussed 
separately  here. 

a.  Large  Impactors 

The  procedures  for  the  clamped  and  simply  supported  plates  are  the 
same  but  different  curves  are  employed.  For  simply  supported  plates,  read 
g^(n)  from  the  lower  curve  of  Fig.  4  and  e  corresponding  to  the  value  of  M 
from  Fig.  26;  for  clamped  plates,  read  g2(n)  from  the  upper  curve  of  Fig.  4 
and  T  from  Fig.  27.  Then,  the  maximum  strain  may  be  computed  from  the  formula 


where 


e  =  g .  (n)  e 

X  ,max  °i 


(105) 


V  =  impact  velocity 
h  =  plate  thickness 
p  =  plate  density 


=  plate  flexural  rigidity 
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Figure  27.  Selected  Design  Curve  for  Clamped 
Plate  Impacted  by  Large  Impactors. 
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10.0  20.0  50.0  100.0 

Mass  Ratio,  M  = 


28.  Selected  Design  Curve  for  Impact 
of  Simply  Supported  Rectangular 
Plate  by  a  Small  Impactor. 
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This  procedure  will  be  Illustrated  by  an  example. 


For  a  rectangular  simply  supported  plate  with  the  following  properties 


Dll  =  680.  N-m 

D22  =  2490.  N-m 
a  =  170  mm 
b  *  350  mm 


h  =  7  mm 
p  =  1.7  g/cm^ 
E  =  8.96  Gpa 
V  “  0.3 


Impacted  by  a  projectile  with  the  following  properties 


R  a  25  mm 
m2  =  1.42  kg 
E  =  2.1x10^  Gpa 


V  =  0.33 

V  =  3.0  m/s 


we  first  confute  the  values  of  the  mass  ratio  and  aspect-orthotropy  ratio 


M  =  0.5 


n  =  0.451 


From  Fig.  4,  we  find  g^^ (0.451)  =  0.26,  and  from  Fig.  26,  e  “  0.91  corresponds 
to  M  =  0.5.  The  maximum  strain  In  the  x-dlrectlon  Is 


e  =  (0.26) (0.91) 

x,max 


(3.0) (0.007) 
680 

(1700) (0.007) 


=  6.57  X  10 


b.  Small  Impactors 

The  design  curve  for  Impact  of  plates  by  small  Impactors  can  handle 
only  simply  supported  boimdary  conditions.  To  use  the  design  curve,  read  the 

ic 

value  of  generalized  strain  e  from  Fig.  28  ,  and  the  value  of  g(8  ,n)  from 

X  p 

Fig.  22.  (Note  that  g(8  ,ti)  does  not  depend  significantly  on  3  for  3  >7). 

p  P  P 

Then,  the  maximum  strain  may  be  coi!9>uted  from  the  formula 


^x,max  =  ^x 


(106) 
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where 


2  pr- 

$  =  0.09897 


11  V  / 


This  procedure  is  illustrated  by  an  example. 

Example  #4 

Consider  the  same  impact  situation  as  in  Example  #3  but  change  the 
impactor  mass  to  m2  =  0.035  kg.  Now  we  have  M=20,  a  small  impactor  case. 
From  Fig.  28,  we  read  e*  =  0.025;  and  from  Fig.  22,  g(e  ,0.451)  =  0.35. 

X  p 

Computing  8p,  we  get 


fO  09897^  (Q-170)  J(1700)  (0.0^ 

(.u.uyay/;  3^+0  453^  W  ^gQ 

^  Q  2  _  1/5 

r (3.0)  (2.008  xio  ) 
(0.035)^ 


=  6.54 


Therefore,  the  maximtim  strain  is 


^x,max  "  ^(0.025) (1.068) (0.35) (6.54) 


(0.007) (3.0) 

f  680 
(1700) (0.007) 


=  6.79  X  10 


B.  Construction  of  New  Design  Curves  by  Experiments 


For  beams  and  plates  with  botmdary  conditions  different  from  those  con¬ 
sidered  in  this  report,  design  curves  may  also  be  constructed  based  on  impact 
experiments.  For  these  structures,  the  parameters  governing  the  impact  re¬ 
sponse  have  the  same  form,  but  the  relationships  (and  thus  the  design  curves) 
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are  different.  For  still  other  structures  (e.g.,  shells,  rods),  new  param¬ 
eters  may  have  to  be  defined. 


1.  Impact  Strain  Curves 

To  construct  new  design  curves  for  large  impactors,  one  may  perform  a 
series  of  liiq>act  experiments  and  measure  the  maximum  bending  strain  which 
occurs  in  the  structure  during  each  Impact.  Then,  compute  the  value  of  the 
generalized  strain  according  to  the  appropriate  formula: 


£  =  e 


max  hv 

c  k 


X  x,inax  hv 


for  beams 


for  plates 


Plot  one  point  for  each  impact  in  e  vs.  M  coordinates.  Logarithmic  scales 
are  most  convenient.  Finally,  estimate  a  curve  to  fit  these  points.  This 
design  curve  may  then  be  used  to  predict  the  response  of  similar  structures 
having  different  dimensions  subject  to  various  impact  conditions,  as  described 
in  Section  A  above. 


2.  Critical  Impact  Velocity  Curves 

In  simple  impact-to-failure  experiments  when  the  strain  is  not  simul¬ 
taneously  measured,  all  that  is  determined  is  the  lowest  impact  velocity  (or 
critical  velocity)  at  which  the  structure  fails.  However,  this  is  sufficient 
to  construct  a  design  curve  useful  in  predicting  failures  due  to  other  impact 
situations. 

For  each  experiment,  the  dimensionless  critical  velocity  v,  defined 
for  plates  in  Eq.  (33)  and  for  beams  in  Eq.  (62),  is  plotted  vs.  the  mass 
ratio  M.  Then  these  points  may  be  connected  to  form  the  design  curve. 
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Note  that  for  different  Impact  conditions,  different  mechanisms  (bending, 
shear,  torsion,  etc.)  may  be  responsible  for  failure^  consequently,  several 
design  curves  may  be  required,  especially  to  completely  describe  a  series  of 
widely  varying  impact  experiments. 

Finally,  note  that  since  the  contact  effects  depend  on  the  impact  velocity, 
through  a  non-linear  relation,  this  approach  is  not  useful  in  the  small  im- 
pactor  domain. 

C.  Construction  of  New  Design  Curves  by  Analytical  Tools 

Using  analytical  methods,  design  curves  for  some  structures  which  are 
not  treated  in  this  report  may  be  constructed.  We  can  offer  a  few  suggestions 
on  approaches  to  follow  in  order  to  construct  such  curves. 

1-  Design  Curves  Based  on  Generalized  One-Degree-of-Freedom  Model 

For  impacts  of  other  structures  in  the  large-lmpactor  regime,  the  gener¬ 
alized  one-degree-of-freedom  model  presented  in  Section  II  may  be  applied. 

To  follow  this  approach,  the  equivalent  structural  stiffness  and  the  strain- 
displacement  factor  d^2  must  be  found;  these  may  be  either  derived  from  the 
exact  solution  of  the  corresponding  static  problem  or  computed  numerically 
using  some  approximate  solution  (such  as  the  finite-difference  or  finite- 
element  methods).  Then  the  equation  of  the  design  curve  may  be  derived 
directly  from  Eq.  (8).  For  beams,  the  form  of  the  generalized  strain  will  be 
the  same  as  in  Eq.  (12);  for  plates,  the  generalized  strain  defined  in  Eq.  (21) 
can  be  employed.  For  other  structures,  a  new  form  may  have  to  be  defined. 
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2.  Design  Curves  Based  on  Impact  Calculations 

Several  dynamic  solution  methods  are  available  for  determining  Impact 
response  of  structures  on  a  problem-by-problem  basis.  These  Include  the 
finite-element  method  (In  particular,  NASTRAN;  see  Ref.  [1]),  Timoshenko- 
type  solution  methods  (application  of  this  approach  to  several  beam  and  plate 
structures  Is  outlined  In  [5]),  and  possibly  some  approximate  solutions  based 
on  the  Timoshenko  method  (similar  to  those  presented  In  this  report  for  small 
Impactors).  Design  curves  may  be  constructed  by  plotting  the  results  of  a 
few  calculations  using  one  of  these  methods  In  coordinates  of  the  appropriate 
generalized  strain  (e  for  large  Impactors,  e*  for  small  Impactors)  vs.  the 
mass  ratio  M  over  the  range  of  M  of  Interest.  These  plotted  points  may  then 
be  connected  to  form  the  design  curve. 
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V.  SUMMARY 

In  this  report,  a  method  for  generating  a  design  curve  which  predicts 
the  peak  response  of  a  structure  subjected  to  low-velocity  impact  has  been 
presented.  Only  a  few  experiments  or  analytical  calculations  are  required 
to  construct  such  a  curve  for  a  given  type  of  structure.  The  importance  of 
such  a  curve  is  that  all  impact  cases,  involving  various  impact  velocities, 
structural  dimensions,  and  material  properties,  fall  on  the  same  curve  within 
a  variation  tolerable  to  the  designer. 

To  facilitate  our  study,  we  have  divided  all  impacts  into  two  regimes, 
large  and  small  impactors,  each  of  which  exhibits  different  kinematic  behavior 
Specifically,  massive  impactors  have  been  observed  to  produce  multiple  impacts 
the  peak  response  of  the  structure  depending  on  the  initial  kinetic  energy  of 
the  impactor;  in  contrast,  small  impactors  strike  the  structure  only  once, 
the  maximum  response  being  governed  chiefly  by  the  impactor ’s  initial  momen¬ 
tum  and  by  contact  effects  at  the  impact  point.  These  differences  in  impact 
response  are  accotinted  for  in  the  distinct  forms  of  the  design  curves  for  each 
domain. 

Design  curves  have  been  developed  here  for  predicting  the  impact  response 
of  simply  supported  beams  and  simply  supported  and  clamped-edge  anisotropic 
plates.  These  curves  have  been  constructed  using  both  results  of  analytical 
solutions  and  data  from  impact  experiments.  Inaddition,  detailed  procedures 
have  been  described  for  using  the  design  curves  presented  in  this  report,  as 
well  as  for  generating  new  curves  for  structures  not  treated  here. 
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APPENDIX  A 


STATIC  SOLUTION  OF  A 
SIMPLY  SUPPORTED  ORTHOTROPIC  PLATE 


In  this  appendix,  we  will  present  the  details  of  the  de¬ 
rivation  of  the  constants  and  ^^2  for  the  central  transverse 
impact  of  a  simply  supported  orthotropic  plate.  As  explained  in 
the  main  text,  these  constants  are  needed  for  developing  the  design 
curve . 


The  static  deflection  of  such  a  plate  due  to  an  arbitrary 
load  distribution  p  (x,y)  is,  according  to  Timoshenko  and  Woinowsky- 
Krieger  M  3]  , 


o» 


a 


9 


(Al) 


where 


and 


For  a  centrally  applied  point  load  P, 


p  a  £n.  <■/„  n2L 

ai  ^  ^  2. 
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The  tensile  strain  in  the  x-*direction  due  to  bending  is 


which  will  have  its  ttiaxitnuin  value  at 


(x,y,z)  =  (a/2,  b/2,  h/2) 


We  Xi7ill  derive  the  relationship  between  the  strain  at  this  point  ^ 
and  the  maximum  deflection. 


e  - 


-  ^ 
^  3  A 


where 


Substituting  for  P,  V7e  obtain 


2 Kin  /£' 

i 


Yf)^/ 


SO  that 


sr 


The  approximation  assumption  used  above,  eq.  (A2),  has  been 
suggested  by  Timoshenko  and  Woinowsky-Krieger  [131  for  simplifying 
the  mathematics  of  several  plate-bending  problems  and  is  exact  for  the 
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case  of  an  isotropic  plate.  The  advantage  of  this  assumption  is  to 
reduce  the  parameters  for  describing  the  geometry  and  anisotropy  of 
the  plate  to  a  single  quantity,  q  •  For  the  composite  plate  on  which 
the  impact  experiments  were  performed,  the  com.puted  values  are 
H  =  Di2  +  2  Dgg  -  1484.6  N~m 

=  1301.2  N-m 

a  difference  of  14%.  « 


/ 
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APPENDIX  B 

STATIC  DEFLECTION  OF  A  CLAMPED  RECTANGULAR 
ORTHOTROPIC  PLATE  DUE  TO  A  CENTPJVL  POINT  LOAD 

In  constructing  an  impact  design  curve  for  a  particular 
structure  according  to  the  generalized  one-degree-of-f reedom  model 
presented  in  the  main  text,  it  is  necessary  to  knox^  the  relation¬ 
ships  between  load,  deflection,  and  strain  in  the  parallel  static 
problem  for  the  structure  of  interest,  Hov7ever,  in  the  case  of  a 
clamped  rectangular  orthotropic  plate,  a  review  of  recent  literature 
indicates  that  an  exact  solution  for  the  static  deflection  due  to 
a  centrally  applied  concentrated  load  does  not  exist.  In  fact,  the 
most  advanced  related  solutions  found  in  the  literature  include  a 
single-term  Ritz  solution  of  the  case  with  uniform  load  (see 
Lekhnitskii  [103  or  Ashton  and  RTiitney  fl43  )  and  an  exact  so¬ 
lution  of  the  simply  supported  case  with  point  load  (see  Am.bartsumyan 
[15]  or  Advanced  Composites  Design  Guide  ii6i  ). 

Therefore,  in  this  appendix,  the  deflection  and  bending 
strain  of  a  centrally  loaded  clamped-edge  rectangular  orthotropic 
plate  are  derived  from  the  solution  for  a  similarly  loaded  and 
supported  isotropic  plate  due  to  Young  [17}  •  This  is  accomplished 

by  reducing  the  governing  equations  for  both  problems  to  the  same 
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dimensionless  form  by  making  a  simplifying  assumption  suggested  by 
Timoshenko  and  Woinowsky-Krieger  Ci33  ,  and  then  transforming  the 
isotropic  solution  into  terms  of  the  dimensionless  variables  corres¬ 
ponding  to  the  orthotropic  case. 


Young’s  Solution 

In  Young’s  solution  for  the  static  deflection  of  a  clamped 
isotropic  plate,  the  solution  for  a  simply  supported  plate  with  a 
central  concentrated  load  is  combined  with  that  for  a  sim.ply  supported 
plate  with  distributed  bending  moments  along  the  edges.  The  edge 
moments  are  then  chosen  so  that  the  deflection  slope  vanishes  at  the 
boundaries.  The  superimposed  solution  consists  of  three  parts, 

w  =  +  W2  +  (bi) 

where  the  first  term  is  the  solution  for  a  simply  supported  rectan¬ 
gular  plate. 


'  -5^  ^  H( 


. . .  '  r'  ' 

^5T 


(B2) 


where 


Z<fL 


The  other  two  terms  are  the  deflections  of  a  simply  supported  plate 
with  moments  applied  along  the  pairs  of  edges  y  =  +  b/2  and  x  =  +  a/2, 
respectively. 


« 


/ 
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27r^D  T, 


mrrx 


•s;.A  ^  ^7 


^  Sn-^p  ^ 


A 

('-/)  S. 


r-'v  ^/n 

— — 77 


r^rrx  .  ^  /pp/r;<  y)  ^  /  yj  /  7 

^“7  - (^ps,  C^sfi  *“J^ 


The  edge  moments  corresponding  to  W2  and  vj^  are 


^rrX 


W,  ...A- r-i)  ^k,  ^ 


in  which  the  coefficients  A^  and  are  determined  from  the  condition 
that  the  slope  at  the  boundaries  is  zero.  Young  has  computed  approxi¬ 
mate  values  of  the  first  few  of  these  coefficients  for  several  values 
of  the  aspect  ratio  (b/a). 
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Modified  Young's  Solution 

In  using  the  solution  for.  static  deflection  to  construct  a 
design  curve  for  impact  response,  we  are  interested  in  the  relation¬ 
ship  between  strain  and  deflection.  In  Young's  solution  as  described 
above,  differentiation  of  the  term  (twice)  to  obtain  an  expression 
for  bending  strain  leads  to  an  infinite  series  which  is  divergent. 

To  avoid  this  difficulty,  an  alternate  solution,  Navier's  solution, 
which  yields  a  bounded  value  of  bending  strain,  is  substituted  for  w^. 
Thus,  we  instead  let 


a 


Transformation  to  Orthotropic  Plate 


The  governing  differential  equation  for  the  deflection  w  of 

V 

an  isotropic  plate  subjected  to  a  lateral  distributed  load  q(x,y)  is 


^4 

d  W 

Ik* 


2 


Letting  x  ^/a,  y  y/b,  and  w  -  w/a,  this  equation  may  be  expressed 


d  w 


'dx 


4/ 
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Likewise,  the  governing  equation  for  an  orthotropic  plate 


//yj 

(BlO) 


may  be  expressed  in  terms  of  the  same  dimensionless  variables  as 


and  where  it  is  assumed  that  H  =  . 

Note  the  similarity  betv7een  eqs.  (B9)  and  (Bll)  •  The  so¬ 
lution  of  Eq.  CB9)  is  of  the  form 


w  =  f(  [a/b]  qa^/D^,  x,y) 


(B13) 


and  the  solution  of  eq.  (B'll)  is  the  same  function  but  involving 
different  dimensionless  parameters, 


w  =  f(  n,  qa^/D^,  x,y) 


(BIA) 


Now  eq.  (B13)  may  be  obtained  by  simply  writing  the  knovvTi  solution, 
eqs.  (Bl),  (B3),  (B4),  and  (b7),  of  the  isotropic  case  in  terms  of 
the  dimensionless  variables  x,  y  and  vj.  Then,  substitution  of  n  for 
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(a/b)^  and  of  qa^/D^  for  qa^/D  yields  eq.  (B14),  which  is  the  solution 
of  eq.  (Bll).  Finally,  by  returning  eq.  (B14)  so  obtained  to  dimen¬ 
sional  variables  x,  y,  and  w,  we  oTstaln  the  solution  of  the  orthotropic 
equation  (BIO) . 

The  terms  to  the  right  of  the  equals  signs  in  eqs.  (B9)  and 
(Bll)  represent  dimensionless  forcing  functions  for  an  arbitrary  dis¬ 
tributed  load  q(x,y),  and  have  the  form  'q  =  qa^/D.  In  the  case  of  a 
uniformly  distributed  load,  q(x,y)  =  q^,  we  would  have  q  =  q^a^/D,  or, 
in  terms  of  total  load  P  =  qosb,  q  =  Pa^/bD.  This  suggests  that  in 
the  case  of  a  concentrated  load  P  the  dimiens ionless  forcing  function  is 
also  q  =  Pa  /bD  for  the  isotropic  case,  or  ^  =  Pa^/bD^  for  the  ortho- 
tropic  case.  Therefore,  in  obtaining  eq.  (B14)  from  eq.  (B13) ,  Pa^/bD 
is  substituted  for  Pa^/bD, 


For  example,  equation  (B7)  may  be  written  in  terms  of  dimen-- 


sionless  variables  as 


(B15) 


This  represents  a  portion  of  the  solution  to  eq.  (B13). 


Substitute  q  for  (a/b)^  and  (Pa2/bD^)  for  (Pa^/bD)  to  obtain 

4F^  ^  ^ /?77rx  <^s/77r7 


Dx  ^ 


^  /?■ 


CB16) 

^  • 

which  represents  a  protion  of  eq.  (B14).  Rewriting  this  equation  in 
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terms  of  dimensional  variables  yields  the  final  solution  for  the  ortho¬ 
tropic  case. 


zr 

^TTx  7)7ry 

i2>S OffS 

(B17) 


Similarly,  from  eqs.  iB3)  and  (BA),  we  obtain 


ly,  - 


rf^rry  .  , 


z: 


&o 


ti9-l 


r /9t7r/  ,  mny  y  ,  / 

ITif  Tif  7^ 


C(^s 

iff 


j99Trk 


(B18) 


J—ys 


MfT} 


By  assembling  eqs.  (B17),  (B18)  and  (F19)  according  to 
eq.  (Bl),  we  obtain  the  solution  for  the  deflection  of  a  rectangular 
clamped-edge  orthotropic  plate  subjected  to  a  centrally  applied  con¬ 
centrated  load.  The  coefficients  and  B^,  which  are  functions  of 
(a/b)  in  the  isotropic  case,  are  now  functions  of 

Evaluating  the  deflection  solution  at  the  center  of  the 
plate  (x  =  o,  y  =  0),  we  obtain  the  maximum  deflection  of  the  plate 


which  may  be  expressed  as 


a 


4/f)  p 


(B20) 
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Combining  eqs.  (B20)  and  (B24)  we  obtain 
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APPENDIX  C 

COMPUTER  PROGRAM  FOR  CALCULATING 
TIMOSHENKO  SOLUTION  OF  PLATE  IMPACT 

This  appendix  presents  a  listing  of  the  FORTRAN  computer  program, 
•named  SMINC3,  for  calculating  the  Timoshenko  small-increment  solution 
of  the  central  transverse  impact  of  a  rectangular  orthotropic  plate, 
Eiq*  ( 40)  • 

The  program  as  listed  is  set  up  to  solve  an  Impact  problem  having 
the  following  parameters: 

V  =  2.45  m/sec 

a  X  b  X  h  =  171  X  349  X  7.1  mm 

m^^  =  0. 660  kg 
m2  =  0.816  kg 

The  plate  bending  stiffnesses,  computed  according  to'  the  Whitney- 
Pagano  (1970)  anisotropic  plate  theory,  are: 

Dll  =  680.4  N-m 
°12  “ 

D22  =  2488.8  N-m 
Dgg  =  505.7  N-m 

°16  “  ®26  “ 

“  2.687  X  10^  N/m 

The  plate  is  treated  as  specially  orthotropic  in  bending  and  the  D 

16 

and  D2^  terms  are  neglected. 

The  program  itself  calculates  the  value  of  the  Hertzian  contact 
stiffness  constant  k2  based  on  supplied  values  of  isotropic  elastic 
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properties  of  the  plate  and  impactor.  The  values  used  are 


plate  (Gr/Ep) 


impactor  (steel) 


Te  =  9.8  X  10^  N/m^ 

=  0.3 


{ 


E  =  2.1  X  10^^ 
V  =  0.33 


N/m^ 


The  contact  radius  of  the  impactor  is  25.4  mm. 
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C  SMINC3  —  PLATE  IMPACT 

COMPUTATION  OF  TIMOSHENKO-TYPE  SMALL-INCREMENT  SOLUTION  OF  HERTZIAN 
IMPACT  OF  A  SIMPLY  SUPPORTED  ORTHOTROPIC  RECTANGULAR  PLATE 
MODIFIED  SUN  €  CHATTOPAOHYAY  METHOD 
DT  -  TIME  INCREMENT 

NDT  -  NO,  OF  TIME  INCREMENTS  CALCULATED 
TOL  -  TOLERANCE  OF  ERROR  IN  NONLINEAR  SOLUTION 
V  -  IMPACT  VELOCITY 

XA*  YB*  H  -  DIMENSIONS  OF  PLATE  IN  X»  Y,  Z  DIRECTIONS 
Ml,  M2  -  MASS  OF  PLATE,  IMPACTOR 
El,  E2  -  ELASTIC  MODULUS  OF  PLATE,  IMPACTOR 
PRl,  PR2  -  POISSON'S  RATIO  OF  PLATE,  IMPACTOR 
NHM,  NHM  -  NUMBER  OF  VIBRATION  MOOES  CONSIDERED 
K2  -  HERTZIAN  CONTACT  STIFFNESS 
XNX,  YNY  -  INITIAL  STRESSES  IN  X,Y  DIRECTIONS 
REAL  L,M1,M2,K2,L11,L12,L13,L22,L23,L33 
DIMENSION  F < 500 0 ),GWI 5000 J ,GXX I  50001 ,GYYI5000) 

PROBLEM  DATA 

NHM=25 

NHN=25 

NDT^lOO 

TOL=1.0E-05 

DT=l,0E-06 

XNX=0. 

YNY=0. 

011=680,4 

012=473.2 

022=2488.8 

066=505.7 

A44=2.687E  07 

A55=A44 

XNX=0, 

yNY=0. 

XA=0.171 

YB=0.349 

H=0.007l 

Ml=0.660 

M2=0.816 

Rl=0.0254 

PR  1=0. 3 

PR2=0.33 

El=9.8E  09 

E2=2.10E  11 

V=2.45 

Z=H/2. 

PI=3. 14159265398 

1001  FORMAT ( 4X, 'N* ,8X, 'TIME* , 14X, 'FORCE* , I 2X, • APPROACH* , I OX, 'DEFLECTION 
1* ,10X, 'STRAIN  X*,10X, 'STRAIN  Y' , lOX, • STRAIN  XY* ) 

DELl={ l.-PRl**2l/El/PI 
0EL2=I l.-PR2**2)/E2/PI 
K2=4.»SQRT(Rl)/(0ELH-0EL2)/PI/3. 

CAY=PI**2/12. 

A44=A44*CAY 
A55=A55*CAY 
B=DT**2/3.0/M2 
C=  4.0/Ml  ,Q2 

P=M1/XA/YB 
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PRINTlOOl 

CALCULATE  NATURAL  FREQUENCIES  OF  PLATE 
D090K=if NOT 
GWIK)=0. 

GXX(K)=0. 

90  GYY(K)=0. 

00247M=1»NHM,2 

TERMA=M*PI /XA 

TERMA2=TERMA1'*2 

L13=A55*T€RMA 

00247N=1»NHN»2 

TERM8=N*PI/YB 

TERMB2=TERMB*»2 

Lll=nil*TERMA2+  D66*TeRMB2*A55 
L12=ID12+D66)*TERMA*TERMB 
L22=066=>TERMA2+  D22>«'TERMB2+A44 
L23=A44*TERMB 

L33=t A55+XNX}*TERMA2+( A44+YNY )*TERMB2 
0=L11*L22-L12**2 

0ETERM=Q*L33+2.*L12*L23*L13-L22*L13**2-L11*L23**2 

AC=I  L12’<'L23-L22*L13)/Q*TERMA 

8C=IL12*L13-L11*L23)/Q*TERMB 

0MEGA2=DETERM/Q/P 

nMEGA=SQRTI0MEGA2) 

IFnM.EQ.l)-AN0-IN.EQ.l))PRINT1002,0MEGA 

1002  FORMAT!*  OMEGAtl,!  =  •tElS.B) 

TERM1=DT*0MEGA 

C2=l.O 

D0247K=1,NDT 

C1=C2 

C2=C0SIK»TERM1) 

AAA=(C1-C2)/0MEGA2 

CONSTRUCT  TABLES  OF  SUMMATION  FUNCTIONS 
GWIK)=GWIK)+AAA 
GXX( K)=GXXCK)+AAA*AC 
247  GYYI  K)=GYYI K)*-AAA*BC 
P=B+C*GW<1) 

Q=l./3. 

S=2,/3. 

D0100N=lfNDT 

CREATE  HEADINGS  FOR  DATA  AT  TOP  OF  EVERY  PAGE 
IFIMOD! Nf61I.NE«0)GOT070 

1003  FORMAT(lHl) 

PRINT1003 

PRINTlOOl 

70  T=N*DT 
A=V»T 

IF{N.EQ.1)G0T014 

BUM=0. 

NM1=N-1 

SUMI=0. 

SUM=0. 

005J=lfNMl 

K=N-J' 

R  =  K 

SUMI=FC J)-SUMI 
SUM=SUM+R*SUMI 
5  8UM=BUM+Ft J)*GWIK+1) 

*1  AO 
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A=A-C*BUM 

A=A-DT**2*l2.*SUM-SUHI/3.}/H2 
COMPUTE  SOLUTION  FOR  FIN)  USING  NEWTON'S  ITERATIVE  METHOD 
GOT013 

14  FNEW=SQRT|A**3) /K2 
13  FN=FNEW 

B=SIGNC l-0tFNI/K2**2-3*0*A*P**2 
FFN=P**3*FN*#3+B*FN**2 
1  '■3.0»A»*2*P^FN-A**3 

FPFN=3«0»P**3*FN*«2<-2.0*B*FN*3.  0<‘A’S'*2*P 
82  FNEW-FN-FFN/FPFN 

ERROR=ABS( IFNEW-FNI/FN) 

IFCERR0R.GT.T0LIG0T013 
IFIN-GT.DGOTOlO 
CHECK  THAT  F(l)  IS  POSITIVE 
IFIFNEW.GT*0«)GOT010 
FNEW=-FNEW 
GOT013 
10  FIN)=FNEW 

IFIFNEW.Lr.0.)F(N)=0. 

CALCULATE  STRAIN  AND  DEFLECTION 
W=0. 

PSIXX=0. 

PSIYY=0* 

D021J=l,N 

K=N-J+1 

W=W«-F(  J)*GWtK) 

PSIXX=PSIXX<-FIJ)*GXXIK) 

21  PSIYY=PSI  YY4-F(  J}*GYYIKi 

w=w«c 

EPSX=-PSIXX*C»Z 
EPSY=-PSI YY»C*Z 
GAMMA=0. 

ALPHA=SIGNI 1»0»FNEW}«( ABS(FNEW}/K2}««S 
PRINTIOOO,N,T,F IN J, ALPHA, W,EPSX,EPSY,GAMMA 
100  CONTINUE 
STOP 

1000  FORMATUX,  14,71  3X,El5,8n 

END 
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APPENDIX  D 

SIMPLY  SUPPORTED  BEAM 
CENTRALLY  LOADED  BY  HALF-SINE  PULSE 


In  this  appendix  are  derived  expressions  for  the  deflection 
and  bending  strain  in  a  simply  supported  beam  subjected  to  a  half¬ 
sine  pulse  laterally  applied  at  midspan. 

The  deflection  of  a  simply  supported  beam  due  to  an  arbitrary 

force  F(t)  applied  at  midspan  is 

ITTX 

(Dl) 


2  A  /^l)  ^ 


J  FM 


where 


2. 


We  will  consider  the  case  where  the  applied  force  is  repre¬ 
sented  by  a  half-sine  pulse  with  respect  to  time. 


F(i) 


7^^  2?czf-<:7r/f 
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Then  the  solution  for  each  time  domain  (during  and  after  the  pulsi 
may  be  written  in  terms  of  the  convolution  integral 

^  ^  ;  iirx 

ft 

f  sin  fr  S/yi  iJi  (t-T)dv,  o<t<  ^-/f 
Jo  ' 

jAjdj  /!  L. 

J  sinfrsm^i/t-rJxfT'^  iQt'  zf->rr/f 

•'o  ■>  '  ^  i 


Evaluating  the  Integral,  the  deflection  during  the  pulse  is  given  by 

w/'x.t)  =  ^  ,(  {!r^ 

frf-  ^  0<Tt<7r/f 

f/^/  (D5) 

If3  =  i  ,  however,  the  denominator  of  this  expression  is  zero,  and 
the  result  is  undefined.  In  this  case,  o)^  =  i^,^.  =  ^ 

steady-state  vibration  problem,  this  would  correspond  to  a  resonant 

forcing  function.  For  this  special  case,  the  following  term  must  be 
substituted. 

‘  l^‘^£X  i-  2.1^  ^ 


—  itos  ft  <ff^sf,ts/ij  2  ft) 
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After  the  pulse,  the  deflection  Is 

7t  S’  (p 

.  ^s/h  Uii-  +  sm  io/t-  Wf)J  ^ 

unless  3  =  i^,  in  which  case  the  ith  term  must  be  replaced  by 

ItL 

^.l^L  sm^  4^  Tr^sft 

'  tt^ET  ^ 


/d!^  &>  57/f 


The  strain  due  to  bending  is 


f=-2 


The  maximum  value  of  this  strain  is  reached  at  midspan.  If  we  define 
a  dimensionless  strain  as 


/\  rr^EI  - 

^  '  h  TL 


(DIO) 


then  the  solution  in  terms  of  this  quantity  is 


J:  4/>»  f/  -  ^  <  ^<rr/f 

'  t  J  (Bii) 


(D12) 
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unless  6  =  1  ,  in  which  case  the  ith  term  must  be  replaced  by 


A 

£  = 


'l  '.3 

jp.(s/n  0<:t<Tr/f 


(D13) 


(D14) 


For  most  problems  involving  small  impactors  (large  M) ,  the  contact 
duration  is  short  compared  with  the  fundamental  period  of  beam  vi¬ 
bration.  Therefore,  it  is  of  interest  to  determine  the  peak  bending 
sbrain  when  the  pulse  duration  is  short.  It  is  observed  in  this  case 
that  the  peak  strain  occurs  after  the  pulse,  that  is,  in  the  domain 
of  eq.  (D12),  which  may  be  written  as 

e  -  z:  fVzJas  ^ 

This  quantity  reaches  its  maximum  at 

2L 


when  it  has  the  value 


p  ^  2. 

Z— *  ^  2/3 

(3 -L 

which  is  a  function  of  B  only.  This  function  is  plotted  in  Fig.  Dl. 

7^£AJ 

UP>  jf-z 
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Figu  Maximum  .  .mensionless  strain  e  vs.  frequency  ratio  6  for  a  simply  supported  beam 

subjected  to  a  half-sine  pulse.' 
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APPENDIX  E 

DETAILED  RESULTS  OF  IMPACT  EXPERIMENTS 


In  this  appendix,  data  from  all  impact  experiments  conducted  at 
Drexel  and  also  some  r^orted  in  the  literature  are  presented  in  both 
tabular  and  graphical  form.  For  canpleteness,  several  figures  which 
appear  in  the  text  have  been  included  in  this  appendix. 


^  •  Results  of  Large-Impactot*  Strain  Measuranent  Experiments 

Strain  measuring  impact  tests  have  been  performed  on  both  aluminum 
and  composite  (  graphite-epojQr  )  plates.  A  total  of  475  experiments 
were  performed  on  ^ecimens  with  clamped  edges  and  5?  3  on  specimen®  with 
sinply  supported  edges. 


Index  to  Large  Impactor  Design  Curves 


Clamped  Plates 


Simply  Supported  Plates 


Figure 

02. 

0.074 

E12 

0.074 

0.111 

E5 

0,111 

0.  l66 

E8 

0,166 

0.189 

El  5 

0.189 

0.227 

E21 

0,227 

0.25 

E3 

0,25. 

0.30 

E18 

0.30 

0.417 

Ell 

0,417 

0,64 

E6 

0.64 

0.80 

El  7 

0.80 

1.0 

El 

1.0 

1.25 

E16 

1.25 

1.56 

E7 

1.56 

2.4 

E10 

2,4 

3.3 

EI9 

3.3 

4.0 

E2 

4.0 

4.4l 

E20 

4.41 

5.29 

El  4 

5.29 

6,02 

E9 

6.02 

9.0 

E4 

9.0 

13.51 

El  3 

13.51 

Figure 

E34 

E27 

E30 

E37 

E43 

E25 

E40 

E33 

E28 

E39 

E23 

E38 

E29 

E32 

E4l 

E24 

E42 

E36 

E31 

E26 

E35 


4 


/ 


110 


TABLE  E1  TRANSVERSE  IMPACT  EXPERIMENTS  ON  CLAMPED  ALUMINUM  PLATES 
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clamped  aluminum  plate  (1/n  ®  0.25). 


Generalized  Strain 


NADC-78259-60- 


Mass  Ratio,  M  =  nij^/m2  ‘ 


—  One-degree-of-freedom  model 
X  Experiment 


Figure  e4#  E'Ciign  curve  and  experimental  data  for  impact  of 
a  clamped  aluminum  plate  (n  =’9.0). 


Generalized  Strain 


N  ADC 78259 -60 


—  One-degree-of-freedom  model 
X  Experiment 


Figure  Design  curve  and  experimental  data  for  impact  of 

clamped  aluminum  plate  (1/n  =  0.111). 


LOGARITHMIC  2x2  CYCLES  46  7200 

KEUFFEL  8t  ESSER  CO.  MADE  in  u.s.A,  ^ 


NADC-78259-60 


NADC-78259.60 


10  1.5  2  2-5  2 _ 4_ 

Figure  E7.  Cyvs.  Mi  Clsuoiqied  Orthotropic  Plate 


6‘  7  8  9  1( 


L35 


LOGARITHMIC  2  x  2  CYCLES  46  7200 

«-  KEUFFEL  &  ES5ER  CO.  made  in  ii  i  a 


LOGARITHMIC  2x2  CYCLES  I  46  7200 

£  KEUFFEL  &  ESSER  CO.  MADF.N  iicA  V 


NADC-7£259-60 


NADC -78259 -60 


Figure  Eli* 


vs*  Ht 


Clamoed  Orthotropic  Plate 


139 


N  ADC -78259  so 


Figure  £12.  Mx  danced  Orthotropxv  rxaze 


140 


UOGARITHMIC  ?.  x  ?.  CYCLES 
KEUFrEL  8c  ESSER  CO  maoc  in  u 


N  ADC -78259 -60 


I 


4 


Figure  £13* 


£yVS.  Mi 


Clashed  Orthotrc^io  Plate 


141 


46  7200 


NADC -78259 -60 


; 

MP 

Figure  E14,^  €^vs.  Ht  Caan^jed  Orthotrcq?ic  Plate 


142 


LOGARITHMIC  2x2  CYCLES 
KEUFFEL  a  ESSER  CO,  MADE  IN  U  S  A. 


N  ADC -78259  <0 


Figure  E15»  ^yvs.  M*  Claimed  Orthotropic  Plate  i 


143 


46  7200 


NADC 78259 -60 


I 

Fij^ure  E17*  i^yvs.  M:  Clan^jed  Orthotropic  Plate 


145 


LOGARITHMIC  2 

KEUFFEL  &  ESSER 


147 


NADC -78259 -60 


5  6  7  3  9  10 


NADC  78259 -60 


rrt- 


TABLE' E3.  ;  TRANSVERSE  IMPACT  EXPERIMENTS  ON  SIMPLY  SUPPORTED  ALUMINUM  PLATES 


NADC -78259 -60 


TABLE  E3»  [continued) 


152 


A 


TABLE  EA 

DATA  FROM  STRAIN-MEASUREMENT  IMPACT  EXPERIMENTS 
ON  SIMPLY  SUPPORTED  GRAPHITE-EPOXY  PLATES 


NADC -78259 -60 


1 

1 

1 

1 

< 

< 

CK 

% 

Os. 

'V 

% 

$ 

♦ 

w 

\S 

to 

VA 

« 

KO 

N 

(0 

CO 

Vn 

Vl 

Ivu 

i« 

>: 

CO 

!? 

N 

N. 

N 

s: 

M 

o 

N 

O 

fo 

'a 

Vo* 

1 

% 

V) 

vS> 

X 

vy 

> 

<  I 

Vi 

(4 

% 

V) 

>• 

rv 

Vo 

0- 

NO 

5$ 

^ 

V 

vn 

>a 

to 

>5 

>8 

Vl 

ro 

Vl 

n 

s 

in 

Di 

o 

(r- 

to 

o 

10 

CO 

lO 

PO 

p> 

? 

O 

‘Sv 

1 

1 

o 

\; 

< 

Vn 

CO 

o 

00 

0^ 

00 

O) 

\S 

0^ 

NS 

M 

to 

>s 

'S  ^ 

K 

§ 

?! 
^  ] 

1 

*M 

% 

s 

Vl 

ts 

<N4i 

i 

v> 

M 

rv 

^1 

o 

sO 

IN. 

<VJ 

a 

C! 

i 

icl 

r« 

CO 

t 

%  1 

1  l!> 

in 

NO 

*n 

1 

VA 

Vl 

*o 

1 

Vo 

1  3  ? 
R  i  i* 

-N, 

< 

<NJ 

0- 

o- 

% 

CM 

•>1 

% 

r>3 

■^1 

2 

*< 

1  1 
1  a  . 

Vl 

CA 

NO 

V\ 

r*i 

NS 

Vl 

% 

N1 

r- 

V 

vl 

1  «  - 
f 

> 

00  1 

r>4  _  ' 

11 

s  ^ 

> 

<3  fl 

153 


TABLE  E4,  continued 


N  ADC  78259 .60 


1 

i 

VO 

• 

o 

II 

cr 

1 
V)  - 

1 

Vb 

ivo* 

o 

>c 

:i5 

o 

>a 

< 

>* 

N# 

n: 

bv 

c>3 

3 

(<? 

^l 

>9 

tsi 

PO 

■.V, 

O'- 

00 

Ivu 

tsl 

< 

rv 

V. 

5? 

< 

PO 

00 

Ov 

% 

N 

o 

O 

% 

oo 

00 

>5^ 

\i 

■ 

v, 

«n 

w 

ea 

^J 

N 

$' 

'a 

N 

•0 

oo 

oo 

v5 

N 

"O 

>a 

X 

s 

M 

!M 

Jo 

1 

§ 

ro 

1 

<o 

\3 

s 

i 

V) 

'\l 

o 

Vo 

s 

V) 

00 

rC 

Ow 

"N. 

!$ 

N 

^sl 

OO 

K 

8 

>; 

8 

V) 

CsJ 

tM 

w 

1 

o 

f^ 

M 

g 

00 

N8 

§ 

Vi 

n 

«o 

0) 

9 1 

1  % 

?>l 

1 

'V. 

V) 

V) 

Vi 

va 

i 

1 

V) 

V) 

?2 

< 

CN 

% 

CN| 

I 

r^ 

'  -c 

V> 

> 

oj 

$: 

•c 

hi 

R  <5  *s 

'O 

va 

Va 

fA 

N3 

r> 

ki 

% 

vs 

fv 

Specimen  B2 

t 

u  5^ 

®0  ‘V. 

^  u 

‘1> 

e  >• 

N 

v5 

u 

154 


i 


TABLE  E4,  continued 


NADC-78259-60 


An  M  ^  tt  >»  ^ 

(0  >»  N  K  (>■ 

^  ^  ^ 

M  CV|  M  fO  (»n  ^ 


V. 


«  s  K 


^  ^  cr>  'h  ^  w  io 

<  c<  ^  ^  to  CO  tn 


I  a  f  S 

p*?  \a  in  ~i  po 


yS 

2 

\o  “ 


8§^s§§§sSSs 

{^SlS^KiSSS  SSs 


w\  >»  in  vn  hft 


Nft  >s  eo  ^  CD  q:>  6> 
V)'oPoc5cJW'^(vi 


>3 

vj  ^ 

X 


8SgSs«8|§|8 

s;?!5SwCS?;m  cmS? 


<V  § 

«U  s' 

f  1-^ 

j  ^ 


\J 

S  ^  '*• 

S  ^5:  -2S. 

^  ^  ? 

Si  > 


^  'O  I 

^  I  ^ 


n|?nj  cm  cm  N  Csl  c4‘'J^j^jjcMj^J 


<C(v}<n'^'>:<nw  - 


<  CM 


V 

hvj  ^  ^ 


155 


TABLE  E4,  continued 


N  ADC  •78259  ^0 


m 

CM 

VO 

• 

o 

II 

cr 


S 

o 

<u 

a 

CO 


I 

'o 

I 


Vb 


! 

I 


’<1 

'■U 

> 

< 


n: 

Vi 


2:^' 

V 


M 

K 


1 1 


I  i  ^ 


52 


w> 

N 

X 


<VJ 


>a 


yi 


tv» 


«V) 


y* 

Nl 


M 


V 

?  r 

I  ^  I 

5  t" 

h  's 


I 


Vo 

NJ> 


(M 

M 


5^ 

V> 


CM 


cvi 


ro 


rr) 

In 


IS. 

Na 


fS! 


tv 


Vo 

rv 

Vi 

M 


0«v 


!c 

rv 


tnI 


OO 


'S. 

#1 

A: 


Vi 

CO 


X 


a 

^1 


Ov 

t>v 


W 


% 


V) 


IV 


«D 


r- 

r- 

cvi 


tv 


Os 

rv 

% 


Vo 

§ 


va 

vi 


eo 

>a 

rt 


o 

o 

<r- 

CO 


ro 


<v 


<SI 


§ 

r- 

M 


i5 


®o 

r» 

CM* 


S: 

<s 

i\ 


V> 

NS 


cvl 


T- 

€Q 


ro 


?> 


N 


^1 

N 

^  Cli 

Vi  II 

.X 


TABLE  E4,  continued 


NADC-78259-60 


157 


NADC-78259-60 


158 


TABLE  E4,  continued 


NAOC -78259 -60 


1 

s 

N 

X. 

Nl 

t 

'S 

0-- 

o* 

% 

Q 

«6 

do 

o 

v» 

C7 

%' 

CO 

5 

$ 

< 

ty 

o 

N 

ri 

« 

IvJ 

Vn 

«o 

o 

o 

d> 

cJ 

xo 

% 

o 

V- 

« 

V. 

"V 

% 

X 

do 

N 

to 

tK 

Si 

1 

Q 

ki 

Of 

S 

5: 

Vj 

H 

CM 

to 

5 

r«^ 

V> 

U) 

M 

3- 

■3“ 

% 

1^ 

0" 

w 

%0 

\5 

2 

\v  ~~ 

X 

1 

8 

CO 

X, 

§ 

I 

X 

1 

N 

N 

■■ 

O 

10 

ar 

a 

80 

N 

s 

?J 

v» 

N 

jr 

> 

d 

0 

Csl 

a- 

Or 

o 

>s 

>» 

V) 

lo 

<0 

s 

JO 

-s 

,X  O 
^  X 

s 

'P 

o 

lO 

r^ 

1 

5i 

1 

§ 

X 

X 

o 

S 

N 

o 

Wl 

% 

fM 

o 

O 

r 

N 

?<  8 
■M  5 

f 

4  % 

*s> 

tsi 

CO 

N 

N 

(M 

(\» 

(M 

(M 

CM 

§  V  ul 

^  ^  t 

^  SS 

> 

% 

rsi 

04 

-Ni 

X 

»o 

3- 

oi 

<T- 

cr 

CNf 

fO 

r- 

iQ 

“X. 

«vi 

. 

‘‘i  ^ 

III 

(^  c 

to 

CO 

\i 

'o 

r* 

IN. 

V 

S  ^ 

^  ^ 

00  ^ 

c>  •> 

vS 

CN 

s  ^ 

'4 

’  ’ 

K  \fs 

^  Vi 

^  4 
_ 

159 


TABLE  E4 ,  continued 


NADC -78259 -60 


1 

8 

N. 

N| 

i 

'S 

t- 

V 

\» 

c\l 

<r 

d 

M 

O 

• 

*o 

n: 

I* 

• 

r 

CO 

o 

N 

cn 

N 

lo 

Q 

N 

IvJ 

rO 

r- 

to 

O 

0- 

£ 

O 

o- 

o- 

«» 

o 

0^ 

5 

fV? 

>c 

to 

«~ 

10 

N 

(M 

rJ 

N 

ri 

1 

Q 

k) 

Of 

5 

5: 

Vj 

>b 

rj 

> 

o 

Vj 

«D 

0“ 

N 

o 

s 

to 

CO 

or 

lO 

o 

1*^ 

2? 

N1 

1 

o 

(« 

SD 

d 

\9 

2 

\u 

X 

10 

tO 

cn 

to 

gr 

8 

N® 

M 

«0 

0 

«0 

w 

o 

M 

Vo 

Si 

M 

>s 

w 

8 

O 

ro 

O 

8 

ro 

in 

r* 

CM 

-3r 

> 

d 

v; 

$ 

ci 

<0 

xS 

r5 

10 

CO 

r- 

04 

o 

r 

'®U 

X 

in 

& 

I 

m 

n 

M 

c 

£! 

1 

8 

k9 

OJ 

o 

8 

N 

1 

rM 

O 

m  5 

J  *»• 

M 

(M 

ro 

Csl 

N 

04 

N 

CO 

(M 

(M 

oJ 

%  V.  UJ 
b  K 

1  S  ^ 
Kf  > 

% 

04 

oi 

?5 

:r 

oi 

CT- 

w 

fO 

r^ 

% 

tv! 

*<  ^ 

^  '0  ? 

1  1  s 
4  ^  c 

U-> 

CO 

'sS 

Vi 

rvj 

'  viS 

^ 

!I^ 

N» 

r>. 

V 

1  $  ' 
N  ^  ^ 

<3  {\ 

JO  ^vv 

1 

K^^HM 

> 


160 


TABLE  E4,  continued 


NADC-78259-60 


4 


Specimen  F2  (r)=0.168) 

( 

1 

5 

'3 

00 

6 

6 

o 

00 

0 

V. 

0^ 

% 

V. 

s 

< 

< 

« 

% 

">s. 

0- 

«0 

« 

IVi? 

1^ 

oo 

o 

o 

<r 

00 

6 

4) 

« 

N 

•>«. 

V 

V. 

V\ 

N 

\ 

eo 

CM 

\ 

cO 

Nik 

c3 

1 

1 

Q 

>u 

5 

5: 

Vi 

5 

> 

o 

o 

O 

CNJ 

r- 

*» 

< 

0 

«o 

% 

V* 

5 

CM 

« 

cr- 

r* 

\9 

\u 

X 

o 

IP 

O 

H? 

m 

>9 

o 

in 

o 

CM 

1 

V! 

«; 

tV) 

o 

ro 

IP 

CM 

Cl* 

CM 

.. 

IP 

CM 

<r^ 

Cm 

o 

lO 

S 

> 

d 

o 

00 

T 

IP 

«o 

>0 

00 

V* 

< 

<>3 

0 

v9 

<0 

6 

.i„^ 

vS 

«0 

O 

0- 

co 

4k 

(M 

.x  a 

Vjj 

m 

M 

€0 

0 

in 

CD 

QO 

2. 

«0 

ro 

:3" 

OS 

In 

i 

£0 

1 

N> 

% 

OJ 

frt 

4) 

s 

GO 

oo 

o 

csl 

8 

S 

uj  R 

f 

CM 

M 

CM 

CNi 

, 

to 

N 

M 

CM 

«Nl 

M 

M 

g  ^  iJ 
$  ^  -2Q. 

:S 

?? 

CM 

<N 

S5 

'Vs 

i  K 

i 

to 
[  ^ 
CM 

t 

cr- 

cr 

w 

. 

. 

I 

?:} 

V 

<Ni 

S  V,  I" 

^  5  § 

f  Q  ^ 

(^  k 

U'v 

CO 

\P 

r^ 

Vi 

'VO 

\D 

r- 

N&‘ 

Tn. 

V 

1  5  ■ 
^  ^ 

CO  ^ 

^  '* 

rr>  ^ 

K  \ 

^  V) 

^  ,1' 

. .  i*i.^ 

161 


TABLE  E4 ,  continued 


NADC-78259-60 


m 

<3* 

<N 

S: 


CO 

a 

0) 

B 


a 

0) 

p. 

CO 


1 

1 

s 

M 

>1 

I 

'S 

iviT 

06 

<3 

tsl 

«0 

% 

o 

!2 

00 

r* 

® 

6 

Os 

N 

tJ 

N 

< 

n 

in 

\ 

N 

S3 

N 

Os 

«M 

1  X 

CM 

in 

00 

o 

r* 

GO 

<0 

>3 

00 

o 

>v 

< 

\3 

< 

Ok 

in 

M 

< 

■S,. 

< 

% 

cn 

0» 

Q 

»u 

>S) 

> 

o 

Vj 

< 

Ci 

. 

o 

1 

0 

o 

1 

(i 

> 

CM 

N 

N 

«o 

V.k 

n9 

Q 

\u  — 

X 

s 

?! 

ca 

^1 

O 

o 

o 

Csl 

\ 

1 

li 

PO 

% 

8 

? 

> 

Q 

>3 

iO 

N 

in 

'0 

"k 

1 

Ci 

> 

>s 

vJ  ^ 

X 

*n 

O 

o 

to 

CO 

;0 

o 

o 

T 

V) 

r^ 

o 

\ 

eo 

1^ 

lo 

o 

M 

1 

1 

«o 

CM 

? 

N 

N» 

1M 

Cd 

M 

M 

CM 

N 

M 

Po 

^  V  ii 

f  ^  ? 

rQ 

> 

CSJ 

4 

04 

»o 

:r 

T- 

vT- 

rv 

rfi 

r^ 

s 

•Vk. 

V 

§  ^  1 

1  Q  ^ 

S  c 

to 

CO 

\i 

r^ 

Vi 

rvy 

M) 

r- 

so* 

In 

1 

V 

s  ^ 

1  3  ' 

- - -  , 

GO 

csi  v 

«VJ  ^ 

^  J‘ 

d  ^ 

•i. 

6x 

Q  U 

I 


162 


TABLE  E4,  continued 


N  ADC -78259 -60 


(K>j^  0^  ^  ^  ^ 

o  O  o  6  ^  U  ^ 


^  £3.  ■j*  o»> 

flQ  ^  ^ 

o  ^  o  o  -C  < 


CM  'W  CK. 

■N,  ^ 


<r-  S\ 


'’  Cn'OjSSSsS^SftSS 


M  — 


#  c 


O  ^ 


51  ^ 


'®o  »  la  £2  ^  S^  J5  ^  ^  § 

^  ^QD^  — 


^  $ 

ki  r' 

f  l*^ 

J  *»• 

^  o 


\J 

g  ^  Ui 

:S  >' 


\>  '0  ^ 

i  I  s 

(5  c 


Ml  M  M 


M  M  (N4  M  M  ^» 


r<\'^  ^  y>  ‘^h^l'o!'^ 

^  K  >  '^.  K  3^1  cr-l  I  is: 


CVJ  W 


ts  ^ 

<  Cm 


I  ^  ' 


163 


TABLE  E4,  continued 


NADC -78259 -60 


4)0  CD  (S  >0 

%  -*  ^  ^  X 

O  ^  ID  'O 


S  »n  ^  ^  _(  > 

*  ^  ^  i>  f!  ’ 


Vo 

> 

tN 

«o  N 
^  o  O  ^ 


O  •?  ^ 

N 


^  ^  5  S  Wi  w  §  > 


<5^  S:r  ^  JSI'  ^ 

r>.  ts^  tV}  V9  Oh  o< 

^  N  M  O  'Q 


>9  «f\«^iO®0«'0®‘DQ^O 

,xo  |iS;^5:^gR§ 


u}  5 

f  l'^ 

^  o 


'J 

A  V  u> 

S  ^  -5G. 

51  ^'' 
>  > 


§  ^  I 
I  ^  ^ 


I  ^  ' 

1^  I"  ^ 


<v\  ^  t:*  ^'>  •'r>  *0  7-  ^ 

rO  I<  >  rv  t<  3-  vT  1* 


<M 


- 1  • 

CJ  W'  ' 


50> 

.r  IN  I  .-j 


I  I  I  s>Lhz  yui  o^t>^  ^  s^'^ 

01^*7  I5*'Z^  O  OOh'Z  O  OQ^  "Z 


TABLE  E4,  continued 


N  ADC -78259 -60 


* 

~ — i 

t 

1 

s 

1 

1 

ivf 

o 

%> 

o 

§ 

> 

V. 

vS 

va 

< 

Os. 

*s. 

> 

'N. 

crs. 

CO 

< 

X 

X. 

.  S*. 

o3 

NO 

«s 

n: 

lN(f 

* 

o 

5: 

«k 

S: 

d 

X.. 

s 

nS 

•S- 

X. 

V) 

Os 

V) 

X 

X. 

•S. 

S 

M 

TABLE  E4,  continued 

Specimen  H2  (ti=0.302) 

1 

1 

• 

% 

v^ 

T 

'•y 

X 

V 

> 

< 

Vi 

> 

% 

>> 

> 

r< 

tK 

Vi 

Vi 

lo 

H 

5n 

ri 

'f  ^ 

xr 

to 

(> 

O 

1 

\ 

In 

tVl 

1 

X. 

s 

(M 

I 

Ox 

X 

\; 

s 

% 

W) 

N 

< 

j 

U) 

< 

fs^ 

Vft 

V 

CO 

0 

wi 

Ci 

Vi 

sa 

i||^HI|S|^s|| 

CO 

o 

to 

5 

00 

p 

Vo 

Os. 

>s 

§ 

>s 

1 

§ 

N 

1 

5 

1 

N 

Vi 

M 

N 

1 1 

1  % 

N 

rsi 

(Vi 

M 

i 

N 

W 

M 

tvr 

N 

111 

“N. 

< 

<N 

>- 

\ 

CM 

> 

♦ 

<N 

■5:: 

csi 

-< 

V 

?  f 
^  ^  1 

1  ^  - 

'o 

CO 

Va 

CO 

Na 

r^ 

5i 

'O 

CO 

% 

va 

fs 

T  ^ 

1  C 

00  ^ 

r>t  vc 

-1 

ts 

?  -  v;s 
yll 
vi 

N  ^ 

165 


TABLE  E4,  continued 


NADC-7€259-60 


66 


TABLE  E4,  continued 


NADC-78259-60 


167 


Generalized  Strain 


NADC -78259 -60 


Mass  Ratio,  M  = 

—  One-degree-of-freedom 

X  Experiment 


i 


Figure 


Design  curve  and  experimental 
supported  plate  (n  =  ^‘O). 


data  for  impact  of  a  simply 
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AT 


0.3  0.5  1.0  2.0 


Mass  Ratio,  M  =  inj^/in2 

-  One— degree— of— freedom  model 
X  Experiment 

Figure  E26*  Design  curve  and  experimental  data  for  impact  of  a 
simply  supported  aluminum  plate  (n  =  9.0). 
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0-5  1.0  2.0 


Mass  Ratio,  M  =  m^/in2 

— .  One-degree-of-freedom  model 

X  Experiment 

Figure  E27#  Design  curve  and  experimental  data  for  impact  of 
simply  supported  aluminum  plate  (l/n  =  0.111). 
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2x2  CYCLES  46  7200 
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Figure  E34#  6^ vs,  Mt  Singly  Supported  Orthotropic  Plato 
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LOGARITHMIC  2x2  CYCLES 
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LOGARITHMIC  2x2  CYCLES  46  7200 

KEUFFEL  8c  ESSER  CO.  MADE  in  u.s  a.  ^ 


loGARITHMIC  2x2  CYCLES  *  46  72^0 
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I  mt 

Figure  E43.  tfy  vs.  Mi  Sin^ply  Supported  Orthotrcraic  Plate 
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Impacted  by  Large  Itapactors. 
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2 •  Results  of  Impact-to-Failure  Experiments 

Each  specimen  was  Impacted  at  a  gradually  Increasing  velocity 
until  damage  was  detected  either  by  inspection  or  by  use  of  a  hand  held 
ultrasonic  transducer  and  oscilloscope  system. 

The  theoretical  failure  strain, used  in  constructing  the  one- 
degree-of- freedom  model  design  curve  is  the  surface  strain  obtained 

Q 

when  the  strain  in  the  outermost  0  lamina  reaches  the  failure  strain 
supplied  by  the  manufacturer  (0.0112). 
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TABLE  'E5- 

IMPACT  FAILURE  EXPERIMENTS  PERFORMED  ON 
COMPOSITE  BEAM  SPECIMENS 


Impact 

Mass, 

©2  (l^g) 

Mass 
Ratio, . 

M=mj^/m2 

Impact 

Velocity, 

V  (m/sec) 

Failure 

Occurred 

Dimensionles! 

Impact 

Velocity, 

V  =  vh/a^ 

A  1  AO  c: 

A  At; A9 

5.63 

No 

0.00333 

o.ioyj 

U«  UDOZ 

6.32 

No 

0.00374 

6.77 

No 

0.00400 

6.95 

No 

0.00411 

7.90 

No 

0.00467 

7.00 

Yes 

0.00414 

7.35 

Yes 

0.00435 

7.53 

Yes 

0.00445 

7.65 

Yes 

0.00452 

8.25 

Yes 

0.00488 

8.66 

Yes 

0.00512 

A  AC^  A 

A  1  AJ^Q 

9.31 

No 

0.00550 

U. UjOD 

• 

U  •  xKJOy 

9.42 

No 

0.00557 

10.30 

No 

0.00609 

11.20 

No 

0.00662 

11.13 

Yes 

0.00658 

11.27 

Yes 

0.00666 

11.47 

Yes 

0.00678 
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TABLE  E5  .  (continued) 


Impact 

Mass, 

in2 

Mass 

Ratio, 

M=m^/ni2 

Impact 

Velocity, 

V  (m/sec) 

Failure 

Occurred 

Dimensionles 

Impact 

Velocity, 

V  =  vh/a^ 

0.0566 

0.1089 

11.72 

Yes 

0.00693 

12.97 

Yes 

0.00767 

0.0280 

0.220 

13.03 

No 

0.00770 

15.45 

No 

0.00913 

17.59 

Yes 

0.01040 

19.17 

Yes 

0.01133 

20.12 

Yes 

0.01190 

23.07 

Yes 

0.01364 

o 

• 

CM 

Yes 

0.01423 

24.10 

Yes 

0.01425 

0.01434 

- 

0.429 

19.05 

No 

0.01126 

21.17 

No 

0.01252 

21.97 

No 

0.01299 

I 

22.66 

Yes 

0.01340 

1 

23.39 

Yes 

0.01383 

26.24 

Yes 

0.01551 

28.42 

Yes 

0.01680 

33.43 

Yes 

0.01976 

192 


Dimensionless  Impact  Velocity, 


NADC-78259-60 


Critical  Velocity  for  Shear  Failure 


“  . Critical  Velocity  for  Bending  Failure 

Experiments  on  Graphite-Epoxy  Beams: 

X  Failure 
•No  Failure 

Figure iE45|#  Critical  impact  velocities  for  shear  and  bending 

failure  and  experimental  data — a  bending-dominated 
case. 
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Figure  E48*  Critical  Impact  Velocity  Curve  for 


Clamped  and  Simply  Supported  Plates 
Impacted  by  Large  Masses. 
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3.  Results  of  Small  Impactor  Strain  Measuring  Experiments 

Plate  specimens  were  Impacted  by  l"(25mm)  radii  aluminum  projectiles 
fired  from  an  air  gun.  Impactor  velocities  were  recorded  using  a  photo¬ 
diode  system.  Plate  strain  was  recorded  using  type  EA-13-062TT  rosettes 
and  a  type  565  dual  beam  oscilloscope. 
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TABLE  E7 

DATA  FROM  EXPERIMENTS  ON  SIMPLY  SUPPORTED 
BEAMS  IMPACTED  BY  SMALL  MASSES 


Beam  Mass, 
Dimensions, 
and  Material 


0.389  kg 

197  X  16.x  16  mm 

steel 


0.601  kg 

305  X  16  X  16  mm 
steel 


0.1096  kg,  aluminum 
152  X  20  X  '12.7  mm 


Mass 

Ratio, 

Impact 

Velocity, 

Contact 

Stiffness, 

M=m^/m2 

V  (m/sec) 

k2(N/ra^^^) 

27.0 

6.37 

21466x10^^ 

It 

13.17 

II 

13.8 

■HIM 

If 

tt 

12.19 

II 

6.85 

5.49 

ft 

II 

8.53 

It 

3.54 

6.28 

ti 

M 

9.33 

If 

41.8 

10.88 

If 

II 

14.87 

11 

23.8 

7.92 

II 

It 

11.06 

It 

11.2 

7.50 

II 

ft 

12.37 

It 

5.88 

6.16 

II 

II 

8.11 

11 

1.93 

4.39 

If 

It 

6.92 

II 

0.0795 


0.0814 


0.1229 


0.1489 


0.2388 


0.2512 


0.3224 


0.3493 


0.0800 


0.0790 


0.1033 


0.1076 


0.1766 


0.1800 


0.2478 


0.2662 


0.475 


0.462 


0.0184 


0.0175 


0.0315 


0.0369 


0.0739 


0.0743 


0.1123 


0.1171 


0.01142 


0.01092 


0. 01703 


0.01716 


0.0341 


0.0330 


0.0554 


0.0579 


0.1689 


0.1569 
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TABLE  E7  (continued) 


Beam  Mass, 
Dimensions, 
and  Material 

Mass 

Ratio, 

M=raj^/m2 

Impact 

Velocity, 

V  (m/sec) 

Contact 

Stiffness 

k2(N/m^^^) 

Generalized 

Strain, 

e  =  ea^/hv 

£ 

0.601  kg 

305  X  16  X  16  mm 

steel 

■ 

41.8 

10.58 

1.744x10^^ 

0.0760 

0.01166 

It 

12.71 

II 

0.0761 

0.01146 

23.8 

9.42 

" 

0.1242 

0.0216 

tf 

12.53 

II 

0.1235 

0.0208 

11.2 

5.33 

II 

0.1534 

0.0328 

tl 

12.01 

11 

0.1611 

0.0317 

5.88 

5.39 

tl 

0.2345 

0.0570 

tl 

7.71 

II 

0.2799 

0.0656 

0.360  kg 

185  X  16  X  16  mm 

steel 

25.0 

11.31 

2.466x10^^ 

0.0868 

0.0203 

tl 

14.42 

II 

0.0978 

0.0223 

12.8 

8.11 

tl 

0.1523 

0.0421 

11 

10.70 

IT 

0.1647  , 

0.0442 

6.34 

8.20 

II 

0.2563 

0.0814 

tt 

10.52 

It 

0.2434 

0.0754 

3.27 

5.73 

II 

0.3252 

0.1221 

tl 

8.60 

II 

0.3203 

0.1155 

0.1074  kg,  steel 

98  X  15  X  9  mm 

1.89 

3.90 

tl 

0.458 

0.2222 

ft 

6.40 

ft 

0.465 

0.2147 
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TABLE  B8  . 


DATA  FROM  GOLDSMITH  (1960) 


Beam  Mass, 
Dimensions, 
and  Material 

Mass 

Ratio, 

M=mj^/m2 

Impact 

Velocity, 

V  (m/sec) 

Contact 

Stiffness: 

k2(N/m^^^) 

Generalized 

Strain, 

e  =  ca^/hv 

c 

3.66  kg 

762  X  25.4  X  25.4 

mm 

steel 

6.44 

0.405 

3.487x10^° 

0.262 

0.0563 

tl 

1.338 

If 

0.286 

0.0500 

It 

2.54 

II 

0.311 

0.0510 

2.12 

0.408 

II 

0.520 

0.1278 

ft 

1.271 

It 

0.564 

0.1237 

If 

—jf 

II 

0.584 

0.1239 

56.9  kg 

2900  X  50  X  50  mm 

steel 

57.0 

2.40 

0.148 

0.00932 

• 

23.7 

2.47 

fl 

0.247 

0.01855 

fl 

3.00 

II 

0*237 

0.01740 

If 

3.40 

If 

0.260 

0.01886 

II 

4.75 

II 

0,239 

0.01676 

11.41 

2.40 

ft 

0.365 

0.03173 

6.34 

2.40 

ff 

0.420 

0.04107 

3.66  kg,  steel 
724x25.4x25.4  mm 

55.0 

1.725 

1.744x10^° 

0.0802 

0.01065 

2.170  kg,  steel 
762x19x19  mm 

266. 

45.7 

1.233x10^° 

0.0268 

0.001486 

1.447  kg,  steel 
762x19x12.7  mm 

177.2 

45.7 

If 

0.0359 

0.965  kg,  steel 
762x12.7x12.7  mm 

118.2 

45.7 

IT 

0.003228 

0.482  kg,  steel 
762x6.3x12.7  mm _ _ 

59.1 

45.7 

-  II  ■' 

' 

0.1578 
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TABLE  E9 

Iiqpact  of  Composite  Plates  by  Small  Masses 


M 

/V“^/FH/S 

P 

/U 

2i.B 

'/.3S 

O,0^7f 

34C0 

0*0^25 

3000 

ZH.& 

3.7<r 

0,0*76’ 

n 

/6ac> 

wam 

s.so 

0,0S6f 

/ioo 

n.i 

3S0 

0,02*77 

^000 

s.8i> 

0,0^^ 

30X> 

2^.7 

3.7S 

0,0€*7 

3000 

2(nB 

2/>s 

OedtB! 

FZ 

32.3 

^000 

^,03 

6, 02*7! 

/lOfi 

18.3 

S.s<f 

6,0/03 

9S0 

18,3 

wm 

EE^I 
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TABLE  ElO 

Data  From  Schwieger 


fiT^  /)^/£ 

;3 

s 

no 

C,2S 

/sn 

0. 072& 

36'^ 

0.50 

Zt.S 

0,03b^ 

53/ 

0,75 

Z33 

0.0368 

730 

/.oo 

2S.0 

C.6377 

757 

A2S 

25.8 

0.033/ 

— 

//SO 

A  SO 

0.0321 

1283 

/.7S 

21. /> 

0.03/3 

WjS 

28.3 

0,0308 

332 

6.SO 

31.8 

0.0237 

HBb 

a.is' 

37.6 

0.0221 

/.OO 

35. /> 

0.02)3 

777 

/.zs 

38.t 

0,0790 

773 

/.so 

3?.7 

0.0/72 

/d/7 

A7S 

0.0 /SI 

nfi 

IQQIII 

0.0  ip 

(  ! 
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4.  Typical  Strain  vs.  Time  Oscilloscope  Traces  for  Impact  of  Clamped  and 
Simply  Supported  Orthotropic  Plates  by  Large  Impactors. 


) 


I 


> 
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a)  Simply  Supported 


b)  Clamped 


Figure  E53:  Impact  of  Specimen  B3,  M  =  0.21,  Impact  Velocity  =  1.73  m/s 
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Scale: 

2000  y  — /div 

wrm 

vertical 

0.5  msec/div 
horizontal 


a)  Simply  Supported  , 

i 


Bi/b/77 


b)  Clamped 


Scale: 

4000  p  ™/div 

iDm 

vertical 
0.5  msec /major 
horizontal 


> 


Figure  E54:  Impact  of  Specimen  B3,  M  =  0.21,  Impact  Velocity  =  2.45  m/sec 
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Scale: 

1000  y  — /dlv 
nua 

vertical 

0.5  msec/div 

horizontal 


a)  Simply  Supported 


I 


Scale: 

1000  V  —/dlv 
mm 

vertical 


0.5  msec/major 
horizontal 


b)  Clamped 


Figure  E55:  Impact  of  Specimen  B7,  M  =  2.38,  Impact  Velocity  =  1.73 
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Scale 

1000  VI  — /div 

Ttmi 

vertical 

0,5  msec/div 
horizontal 


a)  Simply  Supported 


Scale 

1000  p  —/div 

mm 

vertical 


0.5  msec/major 
horizontal 


b)  Clamped 


Figure  E56:  Impact  of  Specimen  B7,  M  =  2.38,  Impact  Velocity  =  2.49  m/sec. 
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Appendix  F 

The  Effect  of  the  Approximation  H  =  on  the 

Predicted  Impact  Response  of  an  Orthotropic  Plate 

To  simplify  the  analysis  of  the  bending  and  vibration  of  specially  ortho¬ 
tropic  plates,  Timoshenko  and  Woinowsky-Krieger  [13]  have  suggested  the  approx¬ 
imation 

®  f  11^^22  (FI) 

where  H  is  defined  as 

»  "  ''12  + 

An  advantage  of  using  this  assumption  is  that  the  equation  governing  bending 
of  an  orthotropic  plate  may  be  transformed,  through  a  simple  change  in  variables, 
into  an  equation  similar  to  the  equation  governing  bending  of  an  isotropic  plate. 
By  reversing  the  transformation,  a  solution  for  the  deflection  of  an  isotropic 
plate  may  be  converted  into  the  corresponding  solution  for  an  orthotropic  plate. 
This  is  done  in  Appendix  B  for  a  centrally  loaded  rectangular  plate  with  clamped 
boundary  conditions;  the  approximation  (FI)  is  Instrumental  in  this  Instance 
because  no  general  solution  exists  for  the  orthotropic  case. 

Like  all  approximations,  however,  application  of  Eq.  (FI)  introduces  some 
error  into  the  solution.  In  this  Appendix,  we  Investigate  the  influence  of  this 
approximation  on  the  strain  predicted  by  the  one-degree-of -freedom  model  of 
impact  of  a  simply  supported  rectangular  orthotropic  plate.  This  is  done  by 
deriving  the  relationship  between  the  generalized  strain  and  the  mass  ratio  M, 
both  with  and  without  the  use  of  Eq.  (FI),  and  then  comparing  the  results. 
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The  static  deflection  of  a  simply  supported  rectangular  orthotropic  plate 

due  to  a  central  concentrated  load  P  is  [10] 

,  mrx  ,  niry  ,  m-rr  .  nir 

sin - sin  — sin  —  sin  — 

a  b  a  b 

■  TT  i  i  T - n — 4 - 

abir  m  n  m_  ^  .  ^..m  n  ,  n_ 

4  11  ^^72  ,4  22 

a  abb 

where  a  and  b  are  the  plate  dimensions  in  the  x  and  y  directions.  If  we  define 
T)  as  in  Eq.  (16)  and 


fll®22 

then  the  maximum  deflection,  which  occurs  at  the  center  of  the  plate,  may  be 


expressed  as 


where 


w  =  P/K, 
max  1 


'  D  b 


Fj^(?,n)  =  I  I  m,n  odd 


m  n  mn 


1.  ^j.o22  42 

B  =m  +2mn  iiC  +  nn 
mn 


Similarly,  the  maximum  x-dlrectlon  bending  strain,  which  also  occurs  at 


the  plate  center,  may  be  expressed  as 


_  2  Pha  _ 

'x,max  2  .  2^  *  ^ 

It  Djjb 


2h  a  K, 


2  „  ■  V* 

’  ‘’ll" 


where 


F,(c,n)  =  I  I 


m,n  odd. 


m  n  mn 
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Therefore,  the  strain-deflection  constant,  defined  in  Eq.  (6),  is 

b  ^ 

*^12  2h  a  F2(C,n) 

To  apply  the  one-degree-of-freedom  inqiact  model  to  this  case,  we  substitute 
Eqs .  (F2)  and  (F3)  Into  Eq,  (8)  with  e  =  0  and  define  the  generalized  strain  e 


as  in  Eq.  (21)  to  obtain 


Gj^(C.n) 


where 


F,(?,n) 

G,(?,ti)  »  ■  ■ 

Note  that  the  approximation  (FI)  corresponds  to  setting  C  =  1.  Thus,  If 
this  assumption  is  used,  then  the  approximate  relation  between  e  and  M. 

X  * 

corresponding  to  Eq.  (22),  is 

G  (l,n) 


Therefore,  the  correction  required  to  adjust  the  approximate  generalized  strain 
given  by  Eq.  (F5)  to  the  exact  value  given  by  Eq.  (F4)  is 

G  (c,n) 

®^^x^  "  G^(l,n)  "  ^ 

This  error  is  plotted  vs.  C  for  a  few  values  of  n  in  Fig.  FI.  Note  that 
the  magnitude  of  this  error  does  not  exceed  17%  for  the  range  of  ?  values 

calculated  and  is  less  than  9%  for  most  lay-ups  commonly  used  in  actual  aircraft 
structures,  which  lie  in  the  range  0.60  <  C  <  1.6,  (Lay-ups  having  ?  values 
outside  this  range  [e.g.,  C  *  0.31  for  unidirectional  AS/3501-6]  rarely  have 
practical  application  as  plate  structures.)  Further,  for  the  range  of  ^  values 
of  the  plate  specimens  used  in  the  impact  experiments,  this  error  is  always  less 
than  6.5%. 
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In  Table  FI,  the  actual  values  of  C,  n,  and  the  error  in  predicted  generalized 
strain  the  plate  specimens  are  suimnarized.  The  largest  value  of  error 

among  all  of  the  plates  is  only  5.4%.  Therefore,  we  conclude  that  the  approxi¬ 
mation  of  Eq.  (FI)  does  not  contribute  a  significant  error  to  the  design  curve 
for  impact  of  simply  supported  plates  by  large  impactors.  We  may  infer  that  a 
similarly  small  correction  is  required  for  orthotropic  plates  with  other  boundary 
conditions . 


Table  FI 

Error  in  Predicted  Generalized  Strain 
for  the  Plate  Impact  Specimens 


Plate 

Number 

5  ^  H 

’^11^22 

.a/  22 

"-<b> 

E(i^) 

B1 

1.430 

0.625 

-3.6% 

B2 

0.164 

-4.2 

B3 

2.39 

-3.7 

B4 

0.625 

-3.6 

B5 

0.074 

-4.0 

B6 

5.29 

-5.4 

B7 

1.38 

-3.2 

B8 

0.282 

-4.1 

B9 

0.625 

-3.6 

FI 

1.252 

0.641 

-2.2 

F2 

0.168 

-2.6 

F3 

2.45 

-2.3 

F4 

0.641 

-2.2 

HI 

1.270 

1.150 

-2.1 

H2 

0.302 

-2.7 

H3 

4.41 

-3.3 

H4 

1.150 

-2.1 
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Error  in  predicted  generalized  strain  vs.  ?  for 
a  few  values  of  the  aspect-orthotropy  ratio  n« 
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